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The thermal decomposition of ethylene oxide exhibits 
an induction period, and is inhibited by the addition of 
inert gases. The assumptions that ethylene oxide mainly 
isomerizes to acetaldehyde, but at the same time produces 
free radicals which induce a chain decomposition of the 


acetaldehyde, and that the recombination of methyl 
radicals requires a third body lead to expressions for the 
length of the induction period and the pressure change 
during a run in qualitative agreement with experiments 
by Heckert and Mack. 





HE thermal decomposition of ethylene oxide, 

studied by Heckert and Mack,' has a 
number of unusual features. The decomposition 
is homogeneous in Pyrex, and follows the 
equation 


(CH2)x0 =CH,+CO. 


There is a pronounced “‘induction period,” lasting 
about thirty minutes at 380°C, which decreases 
as the temperature is increased; the pressure- 
time curve is therefore S shaped. After the 
induction period is over the reaction is approxi- 
mately first order. The addition of inert gases 
increases the induction period considerably, and 
decreases the rate of the reaction as a whole. 
The mechanism proposed by Heckert and Mack, 
the production of an excited acetaldehyde mole- 
cule, which then decomposes to CH,+CO or is 
deactivated by collision, has been shown to be 
improbable by Kassel.” 

Consideration of the structure of ethylene 
oxide indicates that the initial step of Heckert 
and Mack’s mechanism, the conversion of 





' Heckert and Mack, J. Am. Chem. Soc. 51, 2706 (1929). 
_ ? Kassel, Kinetics of Homogeneous Gas Reactions (Chem- 
ical Catalog Co., 1932), p. 284. 
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ethylene oxide to acetaldehyde, is quite plausible, 
especially since an analogous reaction, the con- 
version of cyclopropane to propene, is known.’ 
It is also known that if free radicals are intro- 
duced into acetaldehyde a chain decomposition 
of the acetaldehyde is set up,‘ which proceeds 
rapidly at a temperature considerably lower than 
that at which acetaldehyde will react alone. 
The isomerization of ethylene oxide to acetalde- 
hyde accompanied by the production of a few 
free radicals would, therefore, lead to the ob- 
served reaction. It is the purpose of this paper to 
show that a mechanism of this sort will account 
for many of the features of the decomposition of 
ethylene oxide. 

If we suppose, then, that ethylene oxide 
usually reacts by forming acetaldehyde, but 
that occasionally a free radical is formed which 
reacts with acetaldehyde just as methyls do, we 
can set up the following sequence of reactions: 


3 Chambers and Kistiakowsky, J. Am. Chem. Soc. 56, 
399 (1934). 

4 Allen and Sickman, J. Am. Chem. Soc. 56, 2031 (1934). 

5 All free radicals are denoted by R, and are assumed to 
react in the same way; the main chain carrier presumably 
is, however, the methyl radical produced by II. The chain 
breaking step IV is assumed to require a third body M, 
whose effect is absorbed in the rate constant ky. 
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CH2—CH, 
. + 
O 
CH2—CH, 
- Sw 
O > 
Ill. R+CH;CHO—RH+CH;CO—RH+CO+R_= ky 
IV. 2R+M — X+M ka 


Concentration of (CH2)20 R CH;CHO 
Given by C1 C2 y 


— 2R ky 


— CH;CHO ky’ 


Following Rice and Herzfeld,® a steady state 
equation for the concentration of free radicals 
may be set up, 


des dt =2k\c,—2kyc? = ’ 


C= (Ric1/kRs) i, 


The rate at which ethylene oxide disappears is 
given by 
—dce;/dt=kye,+ky't1= key, (3) 


(4) 


The concentration of ethylene oxide at any 
time ¢ is given by integrating (3), which gives, 
if A is the initial concentration of ethylene 
oxide, 


where 


k=k, +k)’ =k,’ if ki<ky’. 


(5) 


Now acetaldehyde is being formed and decom- 
posed at the same time by reactions II and III; 
its net rate of formation is given by 


dy /dt=ky'¢:—kecxy =k1'c1 — ko( ici /ks) by 


=k,’A exp (—xt)—ke(kiA/k4)! 
Xexp (—xt/2)y, 


C1=A exp (—xt). 


(6) 


by substituting for c, from (5). 

It is now convenient to set exp (—x«t/2)=x, 
which maps the ¢ axis from zero to infinity on 
the x axis from one to zero. Eq. (6) becomes 


dy/dx = —2ki'Ax/x+2key(kiA /R4)*/« (7) 
= —Qx+Py, 


where 


P=2ko(kiA/ka)'/k and Q=2k;'A/x=2A. (8) 


The integral of this equation, giving the con- 
centration of aldehyde at any time, or rather x 
value, during a run, remembering that at ¢=0, 


6F. O. Rice and Herzfeld, J. Am. Chem. Soc. 56, 284 
(1934). 
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x=1 and y=0, is 
y=Q/P?{1+Px—(1+P) exp [—P(1—x) ]}. 


(9) 


The concentration of aldehyde given by this 
equation has a maximum during a run given by 
setting dy/dx=0. The x value at which y reaches 
this maximum is 


x*=1—(1/P) log (1+P). 


Finally, the total pressure p during the run is 
given by 


p=aity+2(A —c1—-y) =2A —-—G1-—-y 
= 2A —Ax*?—2A/P?\1+Px 
—(1+P) exp [—P(1—<) }}. 


This expression is easily seen to have the form 
given by the experimental data, and the latter 
may be used to get an idea of the magnitude of 
the constant P. 

The value of x* is that given by the inflection 
point in the pressure-time curve for the reac- 
tion.” For Heckert and Mack’s' run 168m 
(A =107.42 cm) this is about 30 minutes, and 
with their 2; as x, P=190. 

Furthermore, P should have a positive tem- 
perature coefficient, so that raising the tempera- 
ture will increase P. From (10) we see that in- 
creasing P increases x*, and thus shortens the 
induction period, as is found experimentally. 


(10) 


(11) 


INERT CASES 


So far we have neglected the assumption that 
reaction IV requires a third body and that &, 
depends on the concentration of everything 
present in the reaction vessel, which precludes the 
possibility that (11) will hold exactly. According 
to the scheme presented,’ the only effect of 


7 This is not exactly true, as x* as defined here is the 
inflection point of the p, x curve, which is not identical with 
that of the p, ¢ curve. Strictly, P should be calculated from 
the equation 1—x’ =(1/P) log (1+P+Px’+P?x’)/(2Px’+ 1) ob- 
tained by setting d*p/di? =0, where x’ is the x value of the 
inflection point on the 9, ¢ curve. Since P>>1, 


x!1—(1/P) log (P/2)2x*+(1/P) log 2, 


and P becomes 160 instead of 190 as above. 

8] and II are assumed unimolecular and first order. If 
they have “fallen off” to any considerable extent the 
effect of adding inert gases will be decreased, as I and II 
will then be increased simultaneously with IV. The further 
modifications necessary in this case are easily made. 
Experimentally, P decreases with increasing A but « 
remains almost constant, so that apparently &; has fallen 
off considerably. 
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Fic. 1. Schematic diagram of the concentration of aldehyde 
y asa function of the “time function” x =exp (—x«t/2). 


adding inert gas will be to increase k,, and there- 
fore to decrease P. That decreasing P increases 
the induction period has already been pointed 
out; the general effect of changing P may be seen 
from the following. 

Consider dy/dx for two values of P, P; and P2, 
where P;>P2, all other quantities being the 
same, and let the corresponding values of y be 
yi and ye. Then 


dy,/dx—Pyy.= —Qx, (a) 
dy2/dx —P2y2= —Qx. (b) 


Subtracting (b) from (a) and adding Pyye2 to 
each side we get 


d(yi— ye) /dx —P1i(yi—ye) =(P,—P2)ye. 


Now set z=y1—y2, and multiply each side by 
exp (—P,x). Then 


(dz/dx — Pz) exp (— Px) 
=(P,—P2)yeexp (—Pix). (c) 


The left-hand side of (c) is equal to d[z exp 
(—P,x) ]/dx and the right-hand side is positive 
for O<x<1. The function zexp(—P,x) is 
therefore a monatonic increasing function of x. 
Now when x=1, yi=y2=0, and_ therefore 
z=zexp(—P,\x)=0, so that for all x<1, 
zexp (—P x) is negative, and since exp (— Px) 
is always positive, z= y1—y2<0, or y2>y1. Thus 
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if P is decreased y is increased at all times dur- 
ing the run. Fig. 1 is a schematic diagram of y 
as a function of x for different values of P. An 
increase of y means a corresponding decrease in 
the observed pressure, and adding inert gas 
should therefore decrease the rate of the reac- 
tion, again in agreement with experiment. 


THE INITIAL REACTION 


If, as seems reasonable, we can identity « with 
ky’, we can draw some conclusions as to the 
mechanism of reactions I and II. Heckert and 
Mack? find k;’=10"% exp (—5200/RT)_ sec.—. 
This is definitely in favor of a bond breaking 
mechanism® such as 


CH.—CHse és ha 
~ — —CH2.—CH2.—O-— — CH;CHO; 
O 

the “‘2R’’ of I would then be -CH2.—CH2—O-. 
This mechanism assumes that the ring opens 
only at the C—O bond. An equally plausible 
mechanism would be to assume that the 
—CH.—CH.2—O- radical rearranges at once to 
acetaldehyde, but that occasionally the ring opens 
at the C—C bond, giving —CH2.—O—CH2-. 
This radical might then react as the ‘‘2R”’ of I, 
or it might split off CHO, forming a CHe. In 
any case the mathematical equations would be 
unchanged. 


THE RECOMBINATION OF METHYL RADICALS 


F. O. Rice and Herzfeld® assume that methyls 
recombine by a bimolecular association, while the 
mechanism assumed here requires a three-body 
reaction, which also seems preferable theoreti- 
cally. It should be pointed out, however, that the 
chain decomposition of acetaldehyde induced by 
azomethane‘ seems to be consistent with a bi- 
molecular association. The implications of IV 
requiring a third body in other reactions will be 
discussed elsewhere. 


90. K. Rice and Gershinowitz, J. Chem. Phys. 3, 485 
(1935). 
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Wave-Mechanical Treatment of the LiH Molecule 


Jutian K. Knipp, Research Laboratory of Physics, Harvard University 
(Received January 6, 1936) 


Variational treatments of the LiH molecule and LiH* 
ion have been carried out in which the two electrons of the 
inner shell of the Li atom are represented by Slater wave 
functions and the orbitals of the two valence electrons are 
treated after the method of James and Coolidge (but with- 
out the interelectronic distance as one of the coordinates). 
In this treatment elliptical coordinates are used for the 
expression of the wave function for the outer electrons, 
which takes the form of an exponential times a power 
series. The energy matrix for the LiH molecule contain the 
matrix terms for the LiH* ion, for which calculations were 
also made. The lowest state of the LiH* ion is a *Z state 
with the dissociation products normal hydrogen and ionized 
lithium and is the only one-quantum state. The use of a 
very flexible eleven-term function gives at an internuclear 
distance of three Bohr radii an energy of repulsion of 
0.181 ev. The slope of the potential energy curve at this 
point is found to be —0.791 ev per Bohr radius. The 
lowest state of the ion is surely stable with an equilibrium 


distance somewhat greater than three Bohr radii. The 
energies of the four two-quantum 22 states of LiH* at the 
same internuclear distance are obtained from the next 
four roots of the eleventh-order determinantal equation 
used to obtain the energy of the lowest state. These states 
are found to be quite strongly repulsive at this distance. 
On applying the variational method to the ground state 
of the neutral LiH molecule, it is found that the best 
value for the calculated binding energy that can be ex- 
pected is smaller than the observed value by over a half 
volt. The addition of two terms containing the cosine of 
the difference of the azimuthal angles of the two outer 
electrons improves the calculated energy by 0.13 ev. An 
eight-term wave function has been tried for the *2 state 
of LiH dissociating to the same products as the ground 
state. The energy found is a repulsion of 1.89 ev at a 
distance of three Bohr radii, the equilibrium distance of 
the ground state. 





N this paper! we wish to present variational 
treatments of the LiH molecule and LiH* ion 
based on the method of James and Coolidge. This 
method, first applied* to the ground state of Hg, is 
the most accurate method known for the treat- 
ment of simple molecules. For the three states of 
the H2 molecule to which it has been applied, the 
binding energies obtained are in very satisfactory 
agreement with the experimental values. Hut- 
chisson and Muskat?* have treated the LiH mole- 
cule by wave mechanics using the Heitler-London 
method of atomic wave functions and a simplified 
Hamiltonian, obtaining results in fairly good 
agreement with experiment. However, James‘ 
has shown that the approximations made in such 
a treatment have very little justification and that 
the Heitler-London method, when applied in 
rigorous form, gives binding energies which are 
much too small. 


1 A preliminary report on this work was presented at the 
Pittsburgh meeting of the Am. Phys. Society (Dec. 1934); 
see Phys. Rev. 47, 339A (1935). : 

2 James and Coolidge, Phys. Rev. 43, 588 (1933); J. 
Chem. Phys. 1, 825 (1933); J. Chem. Phys. 3, 129 (1935). 
For _ states see R. Present, J. Chem. Phys. 3, 122 
(1935). 

3 Hutchisson and Muskat, Phys. Rev. 40, 340 (1932). 

4H. M. James, Phys. Rev. 43, 588 (1933); J. Chem. 
Phys. 2, 794 (1934). See also J. H. Van Vleck and A. Sher- 
man, Rev. Mod. Phys. 7, 167 (1935). 


The James-Coolidge method is characterized 
by (1) the use of an extremely flexible wave func- 
tion involving a power series in all the coordinates 
and (2) the recognition of the importance of 
terms containing powers of the interelectronic 
distance. Elliptical coordinates are particularly 
suitable since in He, at least, the density function 
has approximately elliptical form. The terms 
containing the interelectronic distance make the 
function very sensitive to the relative positions of 
the electrons and serve to make it small in those 
parts of configuration space in which the two 
electrons are close together. Such terms were first 
introduced by Hylleraas® in obtaining accurate 
energies for the He atom. In carrying out the 
method a number of terms is taken and the co- 
efficients in the series determined to make the 
total energy a minimum. By the variational 
principle, the energies obtained are upper bounds 
to the true energies of the system. 

The LiH molecule is of particular interest be- 
cause it is one of the simplest of the heteronuclear 
molecules. Two of its electrons are closely bound 
forming a closed shell about the lithium nucleus, 
leaving the other two as valence electrons. The 


5E. A. Hylleraas, Zeits. f. Physik 54, 347 (1929). 
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THEORY OF THE LiH MOLECULE 


unobserved molecule HeH has one less electron, 
but all three of its electrons play the part of 
valence electrons. The question of inner shells has 
already been investigated by James’ in his treat- 
ment of Lis. He finds that they can be represented 
satisfactorily by Slater® one-electron functions. 
With such a treatment only the coordinates of 
the valence electrons appear in the series ex- 
pression. This is a great help in simplifying the 
problem. However because of the mathematical 
complexity which still remains, it is found im- 
practicable to include terms containing the inter- 
electronic distance between the valence electrons. 
Hence, an entirely satisfactory solution is not 
obtained in the cases in which inner shells are 
present. It has been found in the treatment of He 
and Lie that the method without these terms 
gives binding energies which are approximately 
one-half volt too small. Our computations on 
LiH give a difference of the same sign and order 
of magnitude. It is presumed that the difference 
is due to the absence of the above-mentioned 
terms or their equivalent, since such is the case 
with He, where no inner electrons are present. 

Incidental to the solution of the problem for 
the neutral molecule, LiH, is the solution of the 
simpler problem, that of the molecular ion, LiH*. 
This molecule has not been observed experi- 
mentally. A detailed treatment of it is of interest 
because of the accuracy of the results obtainable 
and because these results comprise the only 
information we have at the present time of its 
properties. 


THe Litu1um HyprRIDE ION 


Although, with the exception of H2*, we have 
no experimental evidence of diatomic molecules 
having one-electron bonds,’ we have reason to 
believe from the theory that such molecules exist 
having stable states. That they are not observed 
is due to the fact that they do not have extensive 
discrete band spectra. The molecule H2* has 
been shown to have at least two stable states* and 


* Slater, Phys. Rev. 36, 59 (1930). 

’For a recent discussion of Hz, and H.*+ comparing 
experimental and theoretical results, see O. W. Richardson, 
Proc. Roy. Soc. A152, 503 (1935). 

* Morse and Stueckelberg, Phys. Rev. 33, 932 (1929); 
Teller, Zeits. f. Physik 61, 458 (1930); Hylleraas, Zeits. f. 
Physik 71, 739 (1931). See also discussion by Bethe, 
Handbuch der Physik, Vol. 24 (1933). 
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James’ has found that Li,*+ has a stable ground 
state. 

For purposes of comparison it is instructive to 
consider briefly the states of H2*, the theory of 
which has been worked out by Morse and 
Stueckelberg, Teller and others. The lowest two 
states are *S states and dissociate to a normal 
hydrogen atom and a hydrogen ion. One is 
symmetrical in the nuclei and stable, having a 
binding energy of 2.78 ev and an equilibrium 
distance of 2 Bohr radii, while the other is 
antisymmetrical and is repulsive. There are six 
two-quantum states (dissociating to 2s or 2p 
atoms and an ion) four of which are 2 states. 
One is stable, having a binding energy of 1.35 ev 
and an equilibrium distance of 11-5 Bohr radii. 
One of the two 7II states is also probably stable. 
The behavior of the one-quantum states can be 
understood qualitatively as a resonance between 
two unperturbed atomic states with the accom- 
paniment of a resonance energy which is sub- 
tracted in the case of the symmetrical state and 
added in the case of the antisymmetrical state. 
A similar description of the two-quantum states 
can be made. Here the situation is complicated 
by the fact that there are four interacting states. 
They are successively symmetrical and anti- 
symmentrical and it so happens that only one is 
stable. A more careful picture might include an 
examination of the effect of the three-quantum 
states above. It is assumed that the one- and 
two-quantum states, which are about ten volts 
apart, do not affect each other appreciably. 

LiH* is a heteronuclear molecule but its 
energy levels approach those of H2*+ at infinite 
distances for large quantum numbers. A similar 
behavior can be expected at ordinary distances 
for those high quantum number states since in 
such states the valence electron is so far out that 
the combined field of the lithium nucleus and the 
shell is like that of the hydrogen nucleus. There 
being very little overlapping of the wave func- 
tions, the effect of interchange of valence and 
inner electrons is small. We would thus expect 
every other one of the higher 22 states to take on 
antisymmetrical characteristics, that is to have a 
node near the internuclear line. Since the 
lithium atom has no one-quantum state, there is 


°H. M. James, J. Chem. Phys. 3, 9 (1935), 
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but one one-quantum state of LiH* and it has 
the dissociation products, normal hydrogen and 
ionized lithium. The nearest *2 state is 8.17 ev 
higher at infinite separation and has the dissocia- 
tion products, normal lithium (Li 2s) and ionized 
hydrogen. Pauling’® has pointed out that so great 
a separation implies the absence of large reson- 
ance effects, such as appear in H,*t, and he 
suggests that the lower state will not be stable. 
The four two-quantum states of LiH*, cor- 
responding to a dissociation into Li or H 2s or 2p 
and an ion, are less than 2 ev apart and hence we 
can expect resonance interactions and the 
possibility of stable states. If stable states do 
exist, we would expect their equilibrium distances 
to be large. 

The lithium hydride ion is a system of three 
electrons and a nucleus of charge three, the 
lithium nucleus, separated a distance R from a 
nucleus of charge one, the hydrogen nucleus. 
Two of the electrons are closely bound forming a 
closed shell about the lithium nucleus. The third 
electron travels about one or the other or both of 
the nuclei as the case may be and can be thought 
of as a valence electron. Let distances be ex- 
pressed in terms of the first Bohr radius of 
hydrogen and energies in a unit equal to the 
energy of the ground state of the hydrogen atom. 
The Hamiltonian of the ‘‘clamped”’ molecule in 
these atomic units is the expression 


6 3 6 2 2 
b=—-2(vet-+-)4E—, 
R 1 S; 


1; >i Vij 

where V,? is the Laplacian for electron i, 7; is its 
distance from the lithium nucleus and s; is its 
distance from the hydrogen nucleus. The distance 
from electron 7 to electron 7 is 7;;. 

We suppose that the molecular wave function 
can be built up as the sum of products of one- 
electron functions. These functions describing the 
behavior of single electrons we shall call orbitals. 
If u(1, a) is a valence electron orbital occupied by 
electron 1 with spin a, and g(1, a) and g(1, 8) are 
inner electron orbitals with spin a and 8, respec- 
tively, an arbitrarily normalized wave function 
antisymmetrical in all the electrons is given by 
the determinant 


1 L. Pauling, J. Am. Chem. Soc. 53, 3225 (1931). 


KNIPP 


u(l,a) g(l,a) g(1, 8B) 
a | 2, 2, B)|. (2) 
7 osedvyl’ a) g(2,e) gf " (2, 

u(3,a) g(3,a) g(3, 8) 


For the inner electron orbitals we shall take 
Slater atomic wave functions. The outer electron 
wave function we shall express as an exponential! 
times a power series in elliptical coordinates. All 
azimuthal angles will be omitted from our func- 
tion. This means that not only do we restrict 
ourselves to the consideration of = states, but 
that we also omit a detailed treatment of the 
effects of the interelectronic repulsions on the 
wave function. By proper choice of the energy 
values taken for the separated atoms, some of the 
errors introduced in the energy of the system by 
these approximations will cancel out in comput- 
ing the difference between the molecular energy 
and the total atomic energy. 
Our orbitals.are 


u(1, a) =2Csf(l, a), 


F(1, a) =Ar"ur!e-™15(a, 01), 
g(1,a)=e-'"8(a, 01), g(1, 8B) =e-*"16(B, 1), 


where A= (rits.)/R, Mi (r1 —s,)/R and 6 and k 
are parameters independent of m and j. (a, 1) is 
the spin function and takes the value 1 or 0 as the 
spin coordinate o; of the electron is a or 8. We 
can now write our wave function as a sum of 
determinants, 


v=2LCHy, f(i, a) g(1, a) g(1, B) 
jn 1 fa) aa) a8), 
6S edV)'f(3, @) g(3,a) (3, B) 


Upper bounds to the correct energies are given 
by minimizing the integral T= f'y*hydv subject 
to the condition for normalization 1= /Yy*¥dv. 
The integrals are taken over all configuration 
space and summed over the spins. As usual in this 
type of problem, we obtain a secular equation for 
the determination of the energy values W, 


| (F’|h| F)—W(F'|1| F)| =0, 


where (F’|h| F) and (F’|1| F) are, respectively, 
the matrix elements of h and 1 between the 
functions y; and y,’. Methods for the evaluation 
of these quantities are discussed by James.* 





2 &@® © 6 4 Ob = |+| wh —™|™ 4 = © RD DO =| 4 fF. wet) me © 4 


THEORY OF THE LiH MOLECULE 


Because of the large amount of numerical 
computation involved in obtaining the matrix 
elements for any set of values of the parameters 
R, 6, k, all work was carried out for one set only. 
The internuclear distance was taken as three 
Bohr radii, which is the experimental equilibrium 
distance of the ground state of the neutral mole- 
cule. The parameter k appearing in the Slater 
functions, if determined so as to minimize the 
energy of ionized lithium, must have the value 
2.6875. We have used the value 2.7, which gives 
an energy for Lit differing by 0.00003 atomic 
unit from the energy obtained by using the 
optimum value. This energy differs by 1.5 ev 
from the experimental value because of the very 
simple wave function used. The choice of the 
value of the parameter in the outer electron 
function is not a very critical one since the differ- 
ence between an exponential having a very good 
value of 6 and the one used can in our case be 
expanded in a series and absorbed by the co- 
efficients C; of the general series (4). 6 is taken 
equal to 1.35 since that is the value used in the 
work on LiH. 

We can safely assume, because of the diffuse- 
ness of our outer electron function, that there is 
no real improvement in the wave functions of the 
inner electrons in the molecule over those in a 
similar treatment of the atom or ion and there- 
fore that the error in the energy introduced by 
the use of simple functions for those electrons 
remains essentially the same in all three cases and 
will cancel out. For states having normal lithium 
as one of the dissociation products, it is preferable 
to take the total energy of the atom as calculated 
by Wilson" (—14.8384Rh) using a four-param- 
eter variational method, since the simplifications 
made by him in the atomic problem are of the 
same nature as those made here. We take this 
value instead of the experimental one because in 
calculating the difference of molecular and atomic 
energies, to use an exact value of the atomic 
energy is to include the large total error of the 
computed molecular energy in that difference. 
To obtain atomic energies for the higher states, 
for which computations are lacking, we have 
added to the computed core energy the observed 
ionization energy for the atomic states—a step 


" E. B. Wilson, Jr., J. Chem. Phys. 1, 210 (1933). 
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justified by the considerable accuracy obtainable 
in computations of ionization energies, even when 
poor core functions are employed. 

The matrix elements have been calculated for 
all functions with m=0, 1, 2 and 7=0, 1, 2, 3 
except the one having m=2 and j=3. The 
secular equation is an eleventh order determi- 
nantal equation the roots of which are upper 
bounds to the lowest 2 energies of LiH*. 

Having calculated the energy at a point, we 
can obtain the slope of the energy curve with 
little extra work. Consider R’=Rq where g is 
nearly one, and let k’=k/g and 6’=6. Since the 
energy of the inner shell is a minimum with 
respect to k, a small change in the value of & is 
negligible. Taking the same value of 6 at our new 
point is of no harm since, as we have seen, the 
exact value of 6 is not important. k enters in the 
expressions for our integrals only in the product 
kR, which remains unchanged. All integrals of 
unity contain the factor R*" where n is the num- 
ber of electrons. For all potential energy integrals 
the factor is R*"—', while for all kinetic energy 
integrals it is R*-*. R enters in no other way. 
Since the factor g*" cancels out in obtaining the 
new energy, we retain the old matrix elements of 
unity and in the matrix terms for the energy 
multiply all potential energy integrals by g~! and 
all kinetic energy terms by g~*. From these new 
matrix elements the energy at a distance Rdg 
from the original point can be obtained and used 
to calculate the slope. In the computations g was 
taken equal to 1.01 (dg=0.01). 

From the first root of the secular determinant 
we find that the lowest state of LiH*+ has a re- 
pulsive energy of not over 0.181 ev at a distance 
of three Bohr radii. The slope at this distance is 
—0.791 ev per Bohr radius. Details are given in 
the two tables. Since the energy of repulsion is so 
small and the slope so large, we conclude, in 
contradiction to the prediction of Pauling, that 
the molecule must be stable at larger distances. 
This can hardly be due to inner shell effects, for 
James® has shown that in Li,*+, where the be- 
havior of these shells should be essentially the 
same, they give rise to a repulsion. It would 
therefore appear that in addition to resonance 
there may be polarization effects which play an 
important, or even predominant, part in the 
formation of a stable one-electron bond. 
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An examination of the importance of the 
various terms is of interest. The three terms with 
mj equal to (20), (21), (22) together lower the 
energy by only 0.006 ev showing that the contri- 
butions of terms with higher powers than }? is 
small. The importance of terms with odd powers 
of wis great; (03,13) together lower the energy by 
0.163 ev. However, this is less than one-tenth the 
value of terms with j7=1. We can safely assume 
that terms with 7=5 would not lower the energy 
by more than 0.02 ev. Similar evidence for terms 
with j=4 gives an even smaller amount for such 
terms. Therefore we conclude that 0.16+0.02 ev 
is a close upper bound to the lowest state of LiHt 
at R=3. Wilson’s calculated value for the ioniza- 
tion energy of Li is too small by about 0.05 ev. 
If we suppose this difference is due to the simpli- 
fications made in the treatment of the inner shell 
and make a similar correction for LiHt+, we have 
0.11 ev as an estimated value of the energy of 
repulsion. 

The next four roots give repulsive energies of 
3.885, 4.292, 4.645, and 5.468 ev for the four 
two-quantum states dissociating into Ht+Li 2s, 
H++ Li 2, and H 2s and 2p+Li*, respectively. At 
infinite separation the first is 1.840 ev below the 
second which is 0.143 ev below the other two. At 
R=3 the separations are 2.199, 0.497, and 0.823 
ev respectively, neglecting the fact that some of 
the computed energies may be better than others 
because of better wave functions. The slope for 
the first of these states has been computed and 
found to be large (—2.855 ev per Bohr radius) 
making a stable state at considerably larger 
distances quite probable. 

The first power terms in \ are the most im- 
portant in these two-quantum states. The first 
power in \ corresponds to the first power of 7 
which appears in two-quantum atomic functions. 
An estimate of the decrease in the energy which 
would be caused by more terms would require the 
evaluation of many numerical determinants. This 
has not been done, hence the energies given must 
be regarded as upper bounds. However, because 
even in the case of the fifth root there are essen- 
tially six coefficients at our disposal, we may ex- 
pect the computed energies of the two-quantum 
states, particularly the lower ones, to be rather 


good. 
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Normalized coefficients for a simple eight-term 
function for the lowest state are given in Table I. 
The coefficients for eleven-term functions for all 
five states are also given. These latter coefficients 
have all been calculated from the same de- 
terminantal equation and therefore the cor- 
responding wave functions are strictly orthogonal. 
We define the valence orbitals U=[1/(fg*dV)! ] 
{2C,f+ag}, the constant a being so chosen that 
U is orthogonal to g. If the C; are so chosen that 


TABLE I. Coefficients of normalized wave functions for LiH~. 
Cn; for R=3, 6=1.35, k=2.7. All the states are repulsive 
at this distance. W’ is the energy referred to isolated H* 
and Li* as zero. W,,’ is the value of W’ for R= =. 








LOWEST STATE 


Lowest STATE 


SECOND STATE 





1.97191 
2.22048 
1.17694 
1.73098 


—0.41862 
—0.36791 
—0.23984 
— 0.54866 


— 1.6667 
— 1.6528 


0.187 


W,,', Rh 
W’, Rh 


W’—W,,' ev 


THIRD STATE 


—7.91415 
0.65278 
17.16824 
—0.79811 


8.14604 

— 2.59954 
—20.09876 
— 0.05094 


— 1.94582 
1.39377 
5.86476 


—0.9273 
—0.6102 


4.292 


2.21981 
2.29138 
1.02505 
1.74190 


—0.70165 
—0.45181 
—0.05147 
—0.55681 


0.07181 
0.02338 
—0.05034 


— 1.66667 
— 1.65329 


0.1811 


FourRTH STATE 


—9.15988 
1.17803 
21.73215 
0.62136 


12.76449 
— 2.21312 
— 26.57546 
—0.97776 


—3.61720 
0.66424 
6.67040 


—0.9167 
—0.5736 


4.645 


— 1.78672 
0.11061 
0.72710 
1.44001 


1.90121 
— 0.62244 
— 2.26936 
— 1.21418 


—0.01359 
—0.21317 
0.53653 


— 1.0597 
—0.77271 


3.885 


FIFTH STATE 


—0.21866 
0.77098 
—0.94663 
0.08929 


0.42614 
1.69240 
1.14955 
—0.03391 


—0.63036 
— 1.20923 
1.02214 


—0.9167 
—0.5127 


5.468 








TABLE II. Energies of LiH* for R=3.03. 








LowEstT STATE 


SECOND STATE 





W,,’ 


— 1.66007 


WwW’ — 1.64844 


W’ — W,,’ ev 
A(W’—W,,’)/4R 


0.1574 
—0.791 


— 1.0531 
—0.77244 


3.799 
— 2.855 
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Fic. 1. Valence electron orbitals for the lowest five 22 
states of LiH* plotted along the internuclear line. The 
functions are normalized when multiplied by a factor 
which is about 0.4. The lowest state dissociates into normal 
hydrogen and a lithium ion. The four other states are two- 
quantum states and dissociate into Li 2s+H*, Li 2p+H*, 
and Lit+H 2s and H 24, respectively. Atomic character- 


istics are noticeable in all the states. These functions are 
orthogonal to the inner shell functions and to each other. 











& Wave Function Ure 





Lowest State 














¥ is normalized, then U is normalized. The U’s 
for different states, like the corresponding y’s, are 
orthogonal. In Fig. 1 are graphs showing these 
functions plotted along the internuclear axis. In 
the fifth state an antisymmetrical tendency is in 
evidence. 


THE Lita1umM HypRIDE MOLECULE 


For neutral lithium hydride the Hamiltonian 
of the ‘‘clamped”’ molecule in atomic units is 


We are concerned in this investigation with > 
states only and shall suppose at first that our 
wave functions can be expressed as the sums of 
products of orbitals of the type used in the 
lithium hydride ion. The wave functions, being 
antisymmetrical in all the coordinates as required 
by the Pauli principle, must be symmetrical in 
the space coordinates of the valence electrons for 
singlet states and antisymmetrical in those co- 
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ordinates for triplet states. If f and h are two 
valence electron orbitals, expressions having the 
proper symmetry are the functions 


, = ——————-)_ (-—1)7| f(1, a)h(2, 
Wri as ped ne Lf )h(2, B) 
+h(1, a) f(2, 8) ]g(3, a)g(4, 8B) (6) 


for singlet states and 


1 
4=—————(-1)*P f(1, a) 
. “a seavyir — 


h(2, a)g(3, a)g(4, B) (7) 


for triplet states. The summations are over all 
the permutations of the four electrons and gp is 
the number of inversions of the permutations 
from 1234. This notation is equivalent to the use 
of determinants. We use it here in preference to 
determinants because later we shall introduce a 
function involving the relative coordinates of the 
two outer electrons and this is easily done with 
this system of notation. 

Let us assume complete wave functions of the 


form 
Y=LCrdma, (8) 


Sih 

in which the orbitals are the functions 

f(A, a) = Ai™p12e- 16 (a, 01), 

h(A, a) =A" "e-™16(a, 0:1), (9) 

g(1, a) =e-*"18(a, 03). 
5, y and & are parameters independent of mnjk 
and 6(a, 01) is the spin function with the spin 
coordinate o;. The evaluation of matrix elements 


and adjustment of the parameters Cy, is ac- 
complished by the methods used in the treatment 


. of LiH*. 


If we write the Slater wave functions for the 
H atom and for the outer electron of Li in 
elliptical coordinates and put R=3, we find the 
parameters in the exponentials to be 1.5 for H 
and 1.0 for Li. 1.5 and 1.0 are probably good 
values for our parameters 6 and y. The work of 
Present on excited states of Hz shows that the 
number of terms necessary in the James-Coolidge 
method is much less for states with nonequivalent 
electrons if these parameters are taken not the 
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same. However, because in this case it made 
possible a marked decrease in the necessary 
numerical work, 6 and y were given a common 
value, chosen to be 1.35. It was taken nearer the 
value for hydrogen since the wave function for 
lithium is more naturally built up out of powers 
of \ than is the hydrogen function. As in the 
problem of the ion, the value of k chosen was 2.70. 

For the 'Y ground state the eleven terms 
(mnjk: 1000, 1001, 1002, 1003, 1100, 1101, 1010, 
1011, 0000, 2000, 1111) gave the energy W’ 
= —2.1939 or a binding energy of 1.82 ev 
(W,,’ = — 2.0597) to be compared with the prob- 
able experimental’ value 2.63 ev. The seven 
other terms (1102, 1103, 1020, 1030, 0001, 1012, 
2001) were tried separately and each gave less 
than 0.03 ev to the energy. To estimate the value 
of further terms is hazardous. The changes in the 
energy caused by the terms tried seem to indicate 
that the energy improvement obtainable by 
adding seven or eight more terms to the eleven- 
term function would not be over 0.15 ev. It is 
therefore perhaps safe to say that the limiting 
value for the binding energy obtainable by a 
function with a large number of terms of the 
kind here considered is about 2.0 ev or six or 
seven-tenths of a volt from the experimental 
value. Coolidge and James found that their 
simple power series function without the inter- 
electronic terms came within 0.46 ev of the 
observed value for the ground state of He. 
Present’s calculation, gave a difference of about 
0.37 ev for an excited two-quantum state of He 
using a thirteen-term function. James found that 
for the ground state of Liz the difference of 0.5 
ev remained between the experimental and com- 
puted values for the binding energy after taking 
eighteen terms of the series. 





P(12) = M(12) cos (¢1—¢:), 


M(12) =[(Av?—1) (1 — yi?) (Ao? — 1) (1 — 2?) J}, 
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For the *2 state an eight-term function made 
up of the terms (1000, 1001, 1002, 1003, 1101, 
1010, 1011, 2000) gave an energy W’ = —1.920 or 
a repulsion of 1.89 ev (W.,’=—2.0597) at a 
distance of three Bohr radii. Other terms were 
not tried and the contributions of various terms 
were not examined. Present, working on a two- 
quantum * state of He, found that a six-term 
function of this nature but having two different 
exponential parameters came within 0.08 ev of 
the experimental value. However he also dis- 
covered that when the parameters were given the 
same value the convergence was very poor. The 
orbital antisymmetry of the triplet states, which 
causes the functions to vanish when the valence 
electrons are near each other, explains why 
functions not involving the interelectronic dis- 
tance or the equivalent are better for triplet 
states than for singlet states and hence why the 
interelectronic terms are less important in such 
functions. In our case we can regard the repulsive 
energy 1.89 ev as an upper bound to the energy 
which is probably within a few tenths of a volt of 
the true value. 

It is apparent that further improvements of 
any appreciable amount in the energy for the 
singlet state must come through the introduction 
of terms which account for the relative orienta- 
tion of the valence electrons. We observe that the 
interelectronic distance is expansible in a series 
having in it the powers of the cosine of the 
difference of the azimuthal angles of the two 
electrons. It is proposed next to investigate the 
effect of adding terms containing a cosine factor. 

In order that our function be single valued, we 
must introduce a factor which vanishes properly 
at the poles. Let P(12) denote the function 


(10) 


where q, is the azimuthal angle of electron 1. We can take as additional terms functions of the form 


1 
vara teavine DL, a)h(2, B)+h(1, a) f(2, 8) |P(12)g(3, a)g(4, B). 


The number of observed vibrational levels for the 
ground state of LiH is only five. For this reason an estimate 
of the binding energy involves a long extrapolation. Pro- 
fessor Crawford has kindly informed me that the energy 
we require, using a Morse potential curve and recent 


(11) 


experimental data, is about 2.63+0.05 ev. This is more 
reliable than Nakamura’s earlier value of 2.56 ev. See 
F. H. Crawford and T. Jorgenson, Jr., Phys. Rev. 47, 932 
(1935); G. Nakamura, Zeits. f. Physik 59, 218 (1930). 
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P(12) is symmetric with respect to interchange of electrons 1 and 2; hence the proper symmetry 


conditions are maintained. We define the quantities 


6 2 6 2 6 2 2 
M(1)=[(AP—1)(1— a1?) J, y= —(ve+-+-), Hh=—-2(vet-+-)+—, 
1 r; 


YT, Si 


FA=SH USA), 


Si rie 


F= f—([(fg)/(g?) Je, 


(f’-M\fM) = Sf’ (1)M(1)Hof(1)M(1)d Vit (4/R?) (ff), 


(f'f/PW'h) = S fA) fA) (2/12) P(12)h'(2)h(2)d Vie, 
(f'f/PPH'h) = S f(A) F(1)(2/ri2) P?(12)h'(2)h(2)d Vie, 


(h’hM?) = fh’ (A)h(1) M(1)*d Vi. 


The matrix elements of terms involving P with terms not involving P are 
(F’H’ || fhP) =((F’f/PH'h)+(F’h/PH'f) \(2*), 


We give them for the singlet state only. Matrix elements of P terms with P terms are 


(12) 





(F’H’|1| fhP) =0. 


(f'h’P| | fhP) =(f'h'P|1| fhP)(g*)-* Sg(1)g(2)Aig(1)g(2)d Vie +3 {C(f’ M\fM) (g) 
+2(f’ fM?/g?) —(f'g/Pfg) \(a’hM?)+((f’ M\hAM)(g?) + 2(f'hM?/g) 
—(f’g/Phg) \(h'f M*) +[(h'M\fM)(g*) +2(h’ fM?/g?) — (h'g/P fe) \(f'hM?) (14) 
+[(h’ M\hM)(g?) +2(h'hM?/g*) — (h'g/Phg) \(f’fM*)} 


The terms (1000P, 1001P) were added to the 
eleven-term wave function. The first of the two 
terms gave an improvement of the energy of 0.09 
ev and the second an additional improvement of 
0.04 ev. The thirteen-term wave function gave an 
energy W’ = — 2.2036 or a binding energy of 1.95 
ev. It is hardly possible that further terms of this 
kind would improve the energy by more than 0.1 
ev. The failure of the cosine to give greater im- 
provement probably lies in the fact that it does 
not vary rapidly when the electrons are close 
together.* An exact expansion of the inter- 


‘8 That the wave function must depend almost linearly 
on the interelectronic distance is a result of the work on Ho. 
Consider also the wave function for Ho, yo(1 —beBri2) + - - -, 
Where Yo is the first term in the James-Coolidge function. 
Then if 6 and 6} are estimated by fitting this function to 
their solution, we obtain the approximate values 6=0.16, 
b=0.57. The value of 8 is small showing the effect to be 
a long range one. For the He atom a similar estimate can 
be made using Hylleraas’ wave function. We find 6 =0.18, 


(f'h’P|\1| fhP)=3{ (f' $M") (W’hM?*) + (f’hM?)(h' fM?)} (g?). 


+h (f'f/PPh'h) + (f'h/P?*h'f) }(g*), 








electronic distance involves the use of functions 
with discontinuous derivatives. Higher powers of 
the cosine or terms having cos p(¢:1—¢2), p an 
integer, were not used because their introduction 
would mean an extension of all the tables and 
was regarded as impracticable. Terms having 
cos 2(¢:1—¢2) would probably have greater value 
than those used because the cosine of twice the 
angle varies more rapidly for small angles. 

The writer wishes to express his sincere thanks 
to Professor E. C. Kemble for valuable criticisms 
on this work. He is deeply indebted and very 
grateful to Dr. H. M. James for the benefit of 
many discussions and original suggestions. 


b=0.67. This suggests that, for two electrons which can 
be said to be in the same state except for spin, 6 and b are 
roughly independent of the nuclear system in which the 
electrons are moving. They serve to define a sphere of 
influence of one equivalent electron on the other. 
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Molecular Constants and Potential Energy Curves for Diatomic Molecules. II 


Maurice L. Huaoins, Department of Chemistry, The Johns Hopkins University 
(Received February 19, 1936) 


The potential energy expression U =e~4(*-"12)—(C’e-@"(r-re) previously applied to diatomic 
molecules containing only first row elements is now applied to other types. For most of them 
the same value of a (6.0 10* cm") can be used satisfactorily, the ri2 values then being addi- 
tive and depending only on the row of the periodic table in which the atoms occur. Molecules 
containing but few electrons require different values of a and riz. Values of the equilibrium 
distance 7,, calculated from w,, w-x, and these constants, are in most cases within 0.02 or 
0.03A of the experimentally determined distances (see Table V). The relation between Badger’s 
empirical equation for the calculation of 7, and the equations used in this paper is discussed. 





INTRODUCTION 


N the previous paper of this series! the as- 

sumption that the potential energy of a 

diatomic molecule (not too far from equilibrium) 
obeys an expression of the form 


U(in 10-12 erg) = e@—4(r—112) eo C’e—e' (r—re) (1) 


was applied to the band spectrum data for com- 
pounds containing only elements in the first 
row of the periodic table and containing 8 or 
more electrons. It was found that by taking a 
(in 108 cm) as 6.0 and 72 (in 10-8 cm) as 1.46, 
for all such molecules in all states but highly 
excited ones, values of 7, and of various other 
constants could be calculated which are in quite 
good agreement with the experimental values. 
This agreement is takén to indicate that for 
these molecules the repulsive term in the energy 
expression is practically independent of the 
kinds of atoms in the molecule and of the elec- 
tronic state. The attractive term however 
varies widely from molecule to molecule and from 
state to state. 

An equation, based on the above energy ex- 
pression. was given for the calculation of a 
directly from experimental values of w,, w.x, and 
r,. The calculation is quite sensitive to inaccura- 
cies in the assumed constants however; in 10 
cases out of 36 an imaginary result is obtained. 
The values computed in the other cases average 
close to 6.0, although, as shown in Table I here, 
there seems to be a trend from lower to higher 
values as the number of valence electrons (or 
the sum of the kernel charges) increases. It is 


1 Huggins, J. Chem. Phys. 3, 473 (1935). 


quite possible that the relatively poor agreement 
obtained in calculations of r, for the C2 molecule 
is due to the assumption of too high a value of a. 

In this paper the same type of potential func- 
tion, but, in general, with different values of « 
and fry», is assumed to hold for other diatomic 
molecules than those previously considered. On 
the whole the data are neither so extensive nor so 
accurate as for the molecules containing only 
first row elements, hence the calculations of a 
and fr are less sure. As one consequence, the 
agreement between calculated and experimental 
values of 7,, although on the whole satisfactory, 
is not quite so good as shown in the earlier paper. 


THE CALCULATION OF @, f12 AND f,- 


The results of calculations of rough values of a 
by the method mentioned in the second para- 
graph above are collected in Table I. The equa- 
tion? used is 


a= F/2+(11F?/12—0.4769Mw,x.)}, 
F=aw,/2B2r.+3/r.. 


(2a) 


where (2b) 


M is the reduced mass in atomic weight units. 
The other symbols have their usual significance. 
In the case of the hydrogen molecule 1, values 
were calculated from the ro values by the method 
described in the previous paper, taking a=5.0. 
The assumption of a=4.0 or 6.0 for this purpose 
makes but little difference. 

It is evident from the figures in Table I that, 
if the assumed type of potential function is ap- 


* Eqs. (15) and (26) of the previous article contain 
typographical errors. The correct equations are given here 
as (3a) and (2a). 
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CURVES 


TABLE I.° Calculated values of a. 








SuM OF KERNEL 
CHARGES 


Sum oF KERNEL 
CHARGES 





Average 





SUM OF 

KERNEL 

CHARGES 2 3 4 
NaH1.7 BeH* AlH 3.6 


3. 
3. 
BeH 2. 
3. 
1. 





HO 6.2 


MgH* | 
3. 
MegH 2.: 
CaH 2.5 
2. 


Average tf 2.8 
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5 
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plicable throughout, a gradually increases as 
the strength of bond (in the normal molecule) 
increases. Thus, considering hydrides and other 
molecules separately, a increases as the sum of 
the kernel charges (or the number of valence 
electrons) increases. Also, the hydrides of Zn, 
Cd and Hg (all sub-group elements) seem to have 
a higher value of a than those of Be, Mg and 
Ca, in the corresponding main group of the 
periodic table. There is some indication also of 
a slightly higher value of a for the first row 
elements than for corresponding elements in the 
second and lower rows. 

Using these results as a rough guide, calcula- 
tions of 72 from the experimental r, data have 
been made for various values of a to find the 
value or range of values of a for each molecule 
which leads to nearly the same values of rj for 
all of its different states for which the data are 
available. The equations used are as follows: 


* Most of the data used in the preparation of Tables I, 
III and IV were obtained from Jevons, Report on Band 
Spectra of Diatomic Molecules (University Press, Cam- 
bridge, England, 1932). Results are included ior all 
molecular states for which he gives the requisite data (not 
in parenthesis), except that the 26 states for which the 
data lead to imaginary values of a are omitted from 
Table I. For references, other than to Jevons, for work on 
molecules containing only Be, B, C, N, O and F, see the 
previous paper of this series. In Tables III and IV, although 
the values of r, (calc.)—r. (exp.) are given to 3 decimal 
places (except when less are given for r, (exp.), the last 
figure may be in error by 2 or 3 units, as most of the cal- 
culations were made by slide rule. 


a’ = (2.0625a?+0.7154Mw.x,)' 
.750a, 


C’—C=5.85 X10*Mw,?/aa’, 
C=(C’—C)/(a/a’—1), 
r12—1. = (2.303/a) log C. 


It was found that, as with molecules containing 
only first row elements, the same value of a 
could be used successfully for various different 
molecules. Thus for all hydrides containing from 
6 to 8 valence electrons and other diatomic 
molecules containing from 8 to 14 valence elec- 
trons, the assumption of a=6.0 leads to satis- 
factory agreement between the calculated ry 


TABLE II. Values of a and 12 used in calculating r. 
(Huggins) for Tables III and IV. 








A. Hydrides containing 6 to 8 valence electrons and other molecules 
containing 8 to 14 valence electrons: a =6.0 and riz =n +re 
Element r: (or r2) Element mr (or re) 
H 0.38 X 10-8 cm 
Be, B C N,O,F 0.73 
Al, Si, P, S, Cl 0.96 
Ca, Ti 1.00 
B. Hydrides of Zn, Cd, Hg (also Cu, Ag, Au?): a =4.5 
Molecule (CuH?), ZnH (AgH?), CdH 
riz 1.55 1.66 
C. Hydrides of Be, Mg, Ca (also Li, Na, K, Al?): a =3.5 
Molecule (LiH?), BeH (NaH?), MgH, (AlH?) 
ri2 1.59 1.81 
D. Lis: ¢ =2.0. rpie =2.54 


E. Naz, (K2?): ¢ =1.65. rNa2 =3.01, rK, =3.44 (?) 


(AuH?), HgH 
1.66 


(KH?), CaH 
1.92 


F. He: a =4.5, rag =1.01 or a =5.0, rH, =0.95 or a =6.0, rH, =0.90. 
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TABLE III.? Comparisons of calculated and experimental equilibrium distances. 








FORMULA 


We 
(exp.) 
(cm~!) 


Ve 
(exp.) 


(10-8 cm) 


re(cale.) 
—re(exp.) 

(Huggins) 
(10-8 cm) 


r-(calc.) 
—re(exp.) 

(Badger) 

(10-8 cm) 


REFER- 
ENCE 


FORMULA 


We 
(exp.) 
(cm~!) 


Ye 
(exp.) 
(10-8 cm) 


re(cale.) 
—r,(exp.) 
(Huggins) 
(10-8 cm) 


re(calc.) 
—r-(exp.) 

(Badger) 
(10-8 cm) 


REFER 
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—0.028 
—0.059 
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+0.002 
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—0.033 
—0.045 
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1 
1.009 
0. 
i. 
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1. 
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165.39 
323.86 
250.9 
128.0 
214.26 
212.3 
384.6 
705.8 
718.5 
837.86 
1008.12 
530.6 
451.7 
722.3 
269.69 
253.2 
351.60 
123.79 
115.60 
159.23 
75.00 
92.64 
1476.1 
2220.0 
2087.6 
2058.5 
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1448 
1147.5 
1285 
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values. Moreover, these 72 values were nearly 
the same for practically all compounds satisfying 
these requirements which have the same distribu- 


‘ Kirkpatrick and Salant, Phys. Rev. 48, 945 (1935). 

5 Barwald, L. Herzberg and G. Herzberg, Ann. d. 
Physik 20, 569 (1934). 

6 Crawford and Shurcliff, Phys. Rev. 45, 860 (1934). 

7 Jevons gives ro but not 7,. To obtain the latter, use was 
made of the rough value of a computed by means of Eq. 
(21) of the previous paper. 

§ The ro values given for the SiF molecule by Jevons are 
divided by v2 to accord with the conclusions of Badger 
and Blair, Phys. Rev. 47, 881 (1935). The equilibrium 
distance r, was then calculated from ro as outlined in 
reference 7. 


tion of kernel electrons, e.g., for all molecules 
composed of one first row element and one 
second row element. In addition there is ap- 


9G. Herzberg, Ann. d. Physik (5) 15, 677 (1932) ;_Phys. 


Rev. 40, 313 (1932). 

10 The value of r, used is from Badger, Phys. Rev. 46, 
1025 (1934). Van Dijk and Lameris, Physica 2, 785 (1935), 
obtain 79 =1.73A, which is certainly too low, if the values 
of w,. and w,x, are correct. If the data on P2 and Cl» are 
not in error, calculations of the sort made in this paper 
lead one to expect a value of about 1.90A. Maxwell, 
Hendricks and Mosley, Phys. Rev. 49, 199 (1936), obtain 
1.94+0.03A for the internuclear distance in S, molecules 
by electron diffraction. 

1 Brodersen, Zeits. f. Physik 79, 613 (1932). 
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proximate additivity of these 72 values. From 
the set of ‘“‘basic radii”’ given in Table II one can, 
by addition, calculate r2. values and from these 
and the experimental constants w, and wx, 
compute equilibrium distances r, for a large 
number of molecules. The degree of agreement 
with the experimental equilibrium distances is 
shown in Table III. 

It should be mentioned that some of the data 
on which the calculation of the radii in Table II 
are based are meager and inadequate. In many 
cases also the data are for highly excited states, 
for which, in considering only the first row ele- 
ments, poor agreement was found. Slightly better 
agreement between calculated and experimental 
values of 7, could also be obtained if one gave up 
the assumption of additivity or used different 
values of a for different molecules. For present 
purposes at least, with the limited data available, 
the usefulness of the additivity and constancy of 
@ approximations would seem to make their use 
advisable. 

To obtain reasonable agreement between the 
ne values calculated for different states of other 
hydride molecules and of molecules of the alkali 
metals, smaller values of a have to be assumed. 
The values chosen and the 72 values resulting 
therefrom are listed in Table II. The equilibrium 
distances computed from these are compared 
with the experimentally determined distances 
in Table III. Again the agreement is satisfactory. 

For the hydrogen molecule the results of 
calculations of 7:2 assuming three different values 
of a are given (Table IV). There is really little to 


TABLE IV.2 Comparison of calculated and experimental equi- 
librium distances for H». 








r_(cale.) —r_(exp.) 
(Huggins) 
a=4.5 a=5.0 a=6: 
: we(exp.) a re(exp.) r2=1.01 me=0.95 m2=0.90 
SraTE (em!) (em™) (10-8em) (10cm) (10%em) (10-%em) (10-* em) 


+0.013 +0.003 +0.008 
+.000 +.006 
— .005 
—.001 
— .038 


(Badger) 





lso6pr 311, 2293 
Isc5pr Il, 2309 
Iso4pr 31], 2336.5 
leo3pr 31], 2372.53 
lso3po8E+ = 2.197 
lso2pr 311, 2463 +.06 
Iso2s0 3E5+ 2664.61 1. t j —.001 
(3'M)1Z5 2209 : f J -5 
lso3s01Zo+ 2538 —.08 
lso3d7 IIg 2252.90 —.02 
Iso3do Zot 2455.45 —.07 
2po3po1Z_ 2245.4 —.63 
lso2pr Ml, 2446.5 —.03 
lse2po1E,,+ 1357.266 +.02 
Is? 155+ (4371) +.010 


+0.037 
+.033 
+.026 
+.030 


+.06 
—.031 








* The values of a marked with an asterisk were calculated by means 
of Eq. (24) of the earlier paper in this series. 
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TABLE V. Summary of deviations of calculated from experi- 
mental values of r-. 








+(r.(calc. —r,(exp.)) 


0.025A or 0.026Ato 0.051A to 


less 0.050A 0.150A >0.150A 





Huggins 4 


21 
21 
33 


15 


Badger 63 24 








* All those for which Badger gives the necessary constants. 


choose between them, the general agreement 
being about equally good in each case. Whether 
the large disagreement shown, for any choice of 
a, by some states is due to experimental errors or 
is an indication that the assumed potential 
expression is far from accurate for these cases, 
it does not seem possible at present to decide. 


COMPARISON WITH BADGER’S EQUATION 


Badger” has shown that the equilibrium dis- 
tances for most diatomic molecules in most 
of their electronic states for which the data are 
known may be calculated approximately by the 
simple empirical equation 


re=dij;+C;;3 k., (4) 


in which , is the “bond force constant” and C;; 
and d;; are constants characteristic of the rows of 
the periodic table in which the two elements 
occur. As may be seen from Tables III, IV and V, 
the agreement obtained by using this equation 
with the constants he gives is nearly as good as 
by the method described in this paper. In most 
of the cases of marked disagreement between the 
experimental value of the equilibrium distance 
and that calculated by one of the two methods, 
there is also a similar disagreement when the 
calculation is made by the other method. 

The relation between the two methods can be 
shown in the following way. The force constant 
k., in megadynes per centimeter (the unit used by 
Badger) is related to w,, in cm~, by the equation 


ke = 5.85 X10-*Mw,2. (5) 


This may be combined with Eqs. (3a), (3c) and 
(3d) to give 


2.303 100k, 
re=Nio— log (<—). (6a) 
a a*—aa’ 


2 Badger, J. Chem. Phys. 2, 128 (1934); 3, 710 (1935). 
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which is equivalent to 
2.303 a?—aa’ 
ra [rat log ( ) -K,| 
a 100 
2.303 
+|Ku- 
a 


where K;; is any distance. Now we have shown 
that a is practically a constant (usually 6) for 
diatomic molecules of the same 27 class—com- 
posed of elements from the same two rows of the 
periodic table—provided the sums of the kernel 
charges are not too different. Also (see column 7 
of the table in the previous paper) a’ is almost 
invariably small compared with a; moreover, for 
most states of most molecules of the same 7j type, 
a’ has roughly the same value. Therefore (a? — aa’) 
is the difference between a relatively large num- 
ber (36) and a small, roughly constant number 
(approximately 6 for J, 1 type molecules, zero 
for 1, 2 type molecules, etc.). The quantity in 
the first set of brackets in Eq. (6b) is thus nearly 
constant, for a given type, no matter what value 
K;; has. By choosing a suitable value for K;; the 
quantity (K;;—(2.303/a) log k.) can be made 
closely equal to (C;;'/k.}) over quite a range of 
values of k,. Making such a choice for the range 
of k, values of a given 7j type, it is seen that 
Eq. (6b) is approximately equivalent to Eq. (4). 

The relation between Badger’s constants d;; 
and C;;? is readily obtained from a comparison 
of these equations. Equating corresponding 
portions, 


log | (6b) 


2.303 a? —aa’ 
Ki;=fiet log (~— —d;; (7) 


a 100 


Ci 
and Kij= ia 


k,* a 


2.303 
log ke. (8) 


Eliminating’K;; we obtain 


2.303 a?—aa’ 
di;= Tiot log ( ) 
a 100k, 


Ci;3 


—-—. (9 
“ga ) 
For molecules of first row elements, average 
values of k, and a’ (neglecting the few high 
energy states giving poor agreement between 
experimental equilibrium distances and those 
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calculated by Eq. (3)) are 1.28 and 1.23, re- 
spectively. Substituting these and the values 1.46 
for 712 and 6.0 for a, we calculate 


d,,1=1.214—0.92C%4, 3. (10) 


This equation is nearly satisfied by Badger’s 
values of 0.68 and 0.571 for di, and C},,, re- 
spectively. Similar calculations for other 77 types 
lead to similar results. 


CONCLUSION 


It is evident from the agreement between 
calculated and experimental equilibrium dis- 
tances that the type of potential expression used 
is not a bad approximation for values of r not too 
far from r,. The approximate equivalence of the 
constants a and fy, in the repulsive term for 
different molecules composed of atoms of the 
same two rows of the periodic table (having not 
too different a sum of kernel charges) makes this 
type of potential equation especially useful. 
In this paper application has been made only to 
the calculation of 7, values, but other molecular 
constants can also be roughly calculated by the 
methods described in the previous paper. The 
approximate’ additivity of 712 values is useful in 
extending the calculations to other types of 
molecule. 

Undoubtedly in many of the molecules dealt 
with in this paper the true variation of potential 
energy with distance is more accurately ap- 
proximated by a function of the chosen type 
with somewhat different values of a and ry 
than those used here, but with the present data 
one cannot easily (if at all) determine the best 
values; also the usefulness of the function for 
purposes of prediction is greatest if the same 
values of these constants can be used successfully 
with large groups of molecules and states. 

If the prediction of r, were the only useful 
result of the potential energy expression em- 
ployed in this paper, little advantage could be 
claimed over the empirical equation proposed by 
Badger. From the author’s relation however, 
various other unknown or doubtful molecular 
constants can also be calculated. In addition, 
the whole potential energy curve can be obtained 
with reasonable accuracy from a minimum of 
data. 
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II. Perturbations Due to Nearby Vibrational States 


E. Bricut WILSON, JR.,* Mallinckrodt Chemical Laboratory, Harvard University 
(Received February 24, 1936) 


The use of the Wang secular equation for the interpreta- 
tion of the rotational structure of a vibration-rotation 
band of an asymmetrical top molecule is not justified if 
either of the vibrational states involved in producing the 
band is perturbed by another vibrational state of nearly 
the same energy. The present paper discusses the conditions 
under which such a perturbation is to be expected and 


gives the secular equation which should be used instead 
of the Wang determinant in case the perturbation is 
important. This equation assumes a somewhat simpler 
form when applied to orthorhombic molecules, in which 
case it involves only one parameter in addition to the 
three moments of inertia. 





INTRODUCTION 


N a previous paper! (hereafter called I) the 
Hamiltonian operator for the rotation and 
vibration of a semirigid polyatomic molecule was 
obtained and applied to the problem of the 
energy levels of asymmetrical top molecules. It 
was stated, however, that the perturbation tech- 
nique used is not valid when the vibrational state 
of interest interacts with another vibrational 
state of nearly the same energy. The latter special 
case will be treated in this paper. 


THE HAMILTONIAN OPERATOR 


A lower degree of approximation will be used 
than in I, inasmuch as the phenomenon to be 
treated is more complicated. The results are, 
however, probably sufficiently accurate for 
present practical needs. If in the future it should 
become important to have higher approxima- 
tions, these can be obtained by an extension of 
the methods employed in I. 

The quantities 42g and /, which occur in Eq. 
(33) of I will therefore be replaced by their 
approximate values, with the result that the 
Hamiltonian operator H can be written in the 
form used by several other investigators; 
namely, 


H=}(P.*/Ao+P,?/Bo+P.?/Co) 
= (p:P./A ot pyP, /Bot+p-P ./Co) +8, (1) 


in which P,, P,, and P. are the operators for the 


* Society of Fellows, Harvard University. 

‘ E. B. Wilson, Jr. and J. B. Howard, J. Chem. Phys. 4, 
260 (1936). 

* For a full list of references, see footnote 2 of I. 


components of the total angular momentum 
along the three moving axes x, y, and z fixed to 
the molecule; p., p,, and p, are the operators for 
the corresponding components of vibrational 
angular momentum; Ao, Bo, and Cy are the equi- 
librium moments of inertia; B is the operator 
for the vibrational energy. 


SECULAR EQUATION FOR THE GENERAL CASE 


Part of the infinite secular equation for this 
problem (or also the matrix for H) in terms of the 
functions Wryy is represented in Fig. 1, yz being 
a solution of the wave equation for the sym- 
metrical rigid top while yy is a solution of 
Vy = Wry. ve is a function of the rotational 
coordinates only and involves the quantum 
numbers R, whereas yy is a function of the 
internal coordinates and involves the quantum 
numbers V. 

The doubly shaded blocks, representing groups 
of elements diagonal in V but having all values 
of R, R’, come from the purely rotational and 
purely vibrational parts of H in Eq. (1), while 
the singly shaded and unshaded blocks come 
from the coupling terms —p,P,/A,, etc. These 
coupling terms do not contribute to the elements 
diagonal in V because?* 


SWvpvdr=0, 


for nondegenerate states yy. The coupling terms 
may be regarded as a perturbation whose first- 
order effect vanishes. Ordinarily, the second- 
order effect will be fairly small. If, however, two 
states V and V’ have energies Wy and Wy: 


etc., (2) 


3 For proof see I, footnote 11 and associated text. 
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rather near one another, the second-order effect 
will become large, and in fact, the use of ordinary 
perturbation theory will become impractical. 
This case can be treated by retaining the ele- 
ments connecting the states V and V’ and neg- 
lecting all other nondiagonal terms. Thus in Fig. 
1, the unshaded blocks will be neglected but the 
singly shaded elements, which connect two states 
of nearly the same energy, will be retained. The 
secular equation can then be factored into the 
shaded blocks. The energy levels of the two inter- 
acting vibrational states will be determined by 
the double-size, horizontally shaded part of 
Fig. 1. 
The matrix elements for 


(P2*)r, v= Svr*P Pedr, 
etc., are known‘ to be 
(P.*)s, x, wig, K, M=(Py)s, K, Mia, K,M 
= (h? /82?)[J(J+1)—K?], 


(F)s, x, M; 7, &, u=K?*(h?/41?), 
(3) 


(P.*)7, x, mM: 3, K+2, m= —(P,?)y, - M; J, K+2, M 
= (h?/162?)[ (J=K)(J=K—1) 
x (JEK+1)(JEK+2)}}, 


all other elements vanishing. J, K and M are 
the quantum numbers for the symmetrical top 
(represented by R above). With the use of these 
quantities the secular equation can be obtained. 
It has the following nonzero elements *5 


(K, V|K, V)=K?+ov—o; 
(K, V'|K, V')=K?+oy:—o; 
(K, V|K+2, V)=(K, V’|K+2, V’) 
= 2b{[J?—(K+1)*][(J+1)?—(K+1)*}}!; (4) 
(K, V|K, V’)=—(K, V’|K, V)=iG.K; 
(K, V|K+1, V’)=—(K, V’|K+1, V) 
= (iG.+G,)3[(J#K)(J+K+1)]}. 


‘See, for example, O. Klein, Zeits. f. Physik 58, 730 
(1929). . 

5 As in Wang’s solution of the rigid asymmetrical top 
problem (Phys. Rev. 34, 243 (1929)), which this treatment 
somewhat resembles, the secular equation factors into 
equations for each value of J and M; i.e., H is diagonal in 
J and M. It is not diagonal, however, in K. 


BRIGHT WILSON, 
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Fic. 1. Part of the infinite secular equation for H, 
divided into vibrational blocks. Doubly shaded areas are 
diagonal in V. Singly shaded areas connect two vibrational 
states of nearly the same energy. 


This secular equation has 2(2J+1) rows and 
columns, which are labeled by K and V, the 
former having the values J, J—1, ---, —J+1, 
—J, the latter the values V and V’. Conse- 
quently it is of degree 2(2J+1) in o. The symbols 
have the meaning 


a=(1/2c)(1/Ao+1/Bo), 

b=(1/2c)(1/Ao—1/Bo), 

c=1/Co—3(1/Ao+1/Bo), 

o = (82? /h?c)W-—J(J+1)a, oy =(82?/h?c) Wy, 
G,=(2/cAo) > (mana—mini) In, — 


ikl 


G, = (2, ‘cBy) X(n ain — nix) Ti, 


ikl 


G.=(2/cCo) > dxmia—lim i) Ini, 


ikl 
T= S Wv*Q.(0d 0Q,)Pydr. 


The coefficients 1x, mi, ni define the normal 
coordinates Q, by means of the relation 


3N 
Qi = Dd (linéit-minitnig:), (6) 
i=1 
in which &;, n;, and ¢; are the displacements from 
equilibrium in the x, y, z directions, respectively, 
of the ith atom, multiplied by the square root of 
the mass of that atom. If the harmonic oscillator 
approximation is valid, yy and yy can be 
written down and J;,; evaluated (see Eq. 9). 


SECULAR EQUATION FOR ORTHORHOMBIC 
MOLECULES 


Molecules of orthorhombic symmetry (i.¢., 
point groups V, V;, and Coy), a classification 
which includes most asymmetrical top molecules 





LEVELS OF 


of interest, such as water, formaldehyde and 
ethylene, allow additional simplification. For 
this symmetry, p:, p, and p, each has distinctive 
symmetry properties (irreducible representa- 
tions) which are different from the complete 
symmetry of the molecule. Also, G,, for example, 
vanishes unless yyy has the same symmetry 
properties as p,. Therefore, for a given pair V, 
’’, at most one of G,, G, and G., will be different 
from zero. In case py Wy has a symmetry different 
from that of each of p., p, and p., this type of 
perturbation will vanish completely. 

Furthermore, in each of the three cases 
(G.,#0, G,#0, or G.#0) in which the two levels 
may perturb one another, it is possible to factor 
the secular equation into four equations of lower 
degree. Below are given the transformations 
which bring about the factoring and the nonzero 
elements of the factors themselves for the three 
cases. 


Let Ax=27)(Woxut+Ws_xu)Wv 
Ao=Ysomyv, 
Pe=273(Woxut+Ws_xm)Wv' 
To=YPsomyv’, 
ox =2>(Wsxu—Ws_xu)Ww 
Qx«=27(~oxu—Ws_-xu)W' 


(K>0), 


(K>0), 


(K>0), 
(K>0). 


These functions, together with the specification 
of K as even or odd, will serve to factor each of 
the three cases, different combinations occurring 
in each case. 


The elements out of which the various factors 





Case G, #0, Gy=G.=0 
(a) 
(b) 
(c) se 
(d) ox (K ean) 
Case G,¥0, G.-=G.=0 
(a) Ax (K=0, 2, 4, 
(b) Qx (K=2, 4, ---) 
(c) Tx (K=0, 2, 4, 
(d) ox (K=2,4, ---) 
Case G,¥0, G.=G,=0 
(a) Ax and Q« (K=0, 
(b) Ax and Qrx (K= a 
(c) Tx and ¢x (K=0, 
(d) Tx and ox (K=1, 


POLYATOMIC MOLECULES 


+++) and Tx (K=1, a. 
-++) and Ax (K=1, 3, 
and ox (K=1, 3, 
and Qx (K=1, 3,! 


-++) and Qx (K= 
and Ax (K= 
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of the secular equation are constructed are given 
below. 


(Ax | Ax) =K?+6x,1 20J(J+1)+ov—¢ 
(K =0 included), 
($x | ox) =K?—6x,1 3bJ(J+1)+ov—0(K>0), 
(Az| Ao) = (Ao| Ae) = 2-4b[ (J? —1) (J? +27) J}, 
(Ax | Ax+2) = (x | @x+2) = 30{[J?—(K+1)?] 
*X[(J+1)?—(K+1)*]}? (K+1#1), 
(Ax|Qx)=—(Qx«|Ax)=iG.K (G.0), 
(Ao| Ts) = —(T1| Ao) = 2-4%G,[J(J+1) }! 
(G,40), 
(Ax | x41) = —(T'x | Ax+1) = (¢x | 2x41) 
= — («| ox+1) =}1G,[ (JF K)(J+K+1)]}} 
(G,#0) (K>0), 
(Ao| 21) =(Q1| Ao) = 2-9G,[J(J +1) }! (G,#0), 
(Ax | Qx+1) = (Qx41| Ax) 
= +4G,[(J#K)(JtK+1)}! 
(G,#0) (K>0). 
Elements in which Ix replaces Ax and 9x 
replaces Qx are obtained from the above by 
interchanging V and V’ wherever they appear. 
dx,1 is unity if K=1, zero otherwise. 
The combinations of functions which lead to 


the individual factors are given below. The degree 
of each factor is also given. 


-). Degree: J+1. 
. Degree: J+1. 
-). Degree: J. 
. Degree: J. 


-). Degree: J+1. 
--). Degree: J. 

--+), Degree: J+1. 
. Degree: J. 


-++), Degree: J+1 for J even, J for J odd. 
--+). Degree: J for J even, J+1 for J odd. 
+++). Degree: J+1 for J even, J for J odd. 
-++), Degree: J for J even, J+1 for J odd. 





arr oe 


Pay cee 


Seine 
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The matrix elements between functions listed 
under different letters for each case above are all 
zero; hence each letter represents a factor of the 
secular equation. 


QUALITATIVE CONCLUSIONS 


It can be seen from the mathematical treat- 
ment that certain additional conditions must be 
satisfied before two vibrational states V and V’ 
will perturb one another strongly, even when they 
are close together in energy. The first condition 
is that Yy-Wy have the same symmetry properties 
as one or more of the angular momenta p,, p, 
and p,. A second condition is that J;; shall not 
be too small for some pair of normal coordinates 
ki. This requires, since py is approximately given 
by harmonic oscillator functions, that the two 
states V and V’ differ in two and only two of 
their vibrational quantum numbers x and that 
these two quantum numbers change by one unit 
from V to V’. Since Wy 2Wy,, this practically 
limits the state V’ to one in which »,’=v,+1 
and v,’=v,—1 for some pair of normal modes k 
and /. With this restriction, and assuming har- 
monic oscillator functions, 


Tien = 5 (e+ 1) 01(01°/ ,°) J}. (9) 


Finally, there is the third restriction that the 
coefficient of J; in Gz, Gy, or G, must not be too 
small. This coefficient depends on the nature of 
the normal modes of vibration k and / so that it 
is both a geometrical and a dynamical question 
whether a given pair will have a large coefficient 
or not. 

From this discussion it is evident that before 
attempting to analyze the rotational fine struc- 
ture of a given vibration-rotation band of an 
asymmetrical top molecule by the use of the 
Wang equation, it is important to investigate the 
possibility that another vibrational state (which 
may not give an infrared band) perturbs the 
upper state of the band in question. If the quali- 
tative tests given above are satisfied, then it may 
be necessary to use the extended equations given 
in this paper instead of the Wang equation. This 
requires either a knowledge of the exact normal 
modes of vibration or the use of one or more 
additional empirical parameters (G,, etc.). 

Exactly the same situation may arise in sym- 
metrical top molecules if two nondegenerate 
states approximately coincide. 


> 
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A new value for the oxygen parameter in the structure 
of silver phosphate has been determined which gives the 
P—O separation in the phosphate group as 1.61A+0.03A. 
The value found here is slightly larger than that deter- 
mined from an analysis of the crystal KH2PO,. The slight 
increase in P—O distance is explained on the basis of the 
nature of the phosphate group and of some covalent bond 
formation between oxygen and the surrounding silver 


INTRODUCTION 


STUDY of the structure of silver phosphate 
was undertaken to investigate the possible 
influence of bond formation between the silver 
and oxygen atoms on the dimensions of the 
phosphate group. In previous work on this 


* National Research Fellow in Chemistry. 


atoms. To get qualitative agreement between calculated 
and observed intensities it was found necessary to treat 
the silver atoms as oscillators with tetragonal symmetry, 
the symmetry of their environments in the crystal. The 
ratio of the amplitudes of vibration along the tetragonal 
axis and perpendicular to it has been evaluated from the 
data. This result is discussed in relation to the structure of 
silver iodide. 


crystal Wyckoff! has established it as represent- 
ing a cubic type structure shown in Fig. 1. The 
phosphorus and silver atoms are in fixed posi- 
tions and the oxygen positions on the threefold 
axes are determined by one parameter. With 
the methods available at the time Wyckoff was 


1R. W. G. Wyckoff, Am. J. Sci. 10, 107 (1925) ; Zeits. f. 
Krist. 62, 529 (1925). 
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able to make a rough estimate of the oxygen 
parameter, placing it in the neighborhood of 
0.25. This value leads to a phosphorus-to-oxygen 
distance of approximately 2.6A, or nearly 1A 
larger than the same separation observed by J. 
West? in his careful work on KH2PO,. In view 
of the improved methods now existing for the 
interpretation of x-ray intensities it has been 
thought worth while to reinvestigate the struc- 
ture of Ag;PQO, to establish more accurately the 
distortion of the phosphate group suggested by 
the earlier work. 

In the present work a new value for the oxygen 
parameter has been obtained which eliminates 
the large discrepancy mentioned above, and 
indicates only a slight increase in the P—O dis- 
tance over the value found for KH2PO,. In 
addition the interpretation of the x-ray data has 
shown that the silver atoms must be treated 
as oscillators with tetragonal symmetry the 
symmetry of the environment of the atoms in the 
crystal as shown by the dotted lines in Fig. 1. 
For this case a new temperature factor has been 
calculated involving, as an additional parameter, 
the ratio of the force constants along the tetrag- 
onal axis and perpendicular to it. In this way it 
has been possible to account for the data and to 
obtain a rough idea of the potential function of 
the silver atoms in the crystal. This result is 
discussed also in relation to the peculiarities in 
the intensities of reflection from silver iodide* 
in an attempt to give a more satisfactory ex- 
planation of the observed effects. 


? J. West, Zeits. f. Krist. 74, 306 (1930). 
*L. Helmholz, J. Chem. Phys. 3, 740 (1935). 


PROCEDURE 


Single crystals of silver phosphate were ob- 
tained by redissolving precipitated Ag;PQO, in 
concentrated ammonia solution and allowing 
the ammonia to escape slowly. This procedure 
yields small crystals of the yellow solid, the 
largest obtained being about 1.51.50.5 mm. 
The color and habit of these crystals are the 
same as those given by Groth for Ag3;PO, and the 
unit dimensions checked those given by Wyckoff. 
These criteria of chemical composition were 
assumed to be sufficient and the crystals were 
not analyzed. 

Laue and oscillation pictures were taken of 
three different samples from three different 
preparations. The Laue photographs were taken 
by using continuous radiation from a tungsten 
target with Amin=0.24A. For the oscillation 
pictures Mo Ka was used. The relative intensi- 
ties of reflection were obtained by comparing 
lines of equal intensity on photographs of 
different exposure times according to the usual 
practice. A set of five different exposures of 12, 
24, 48, 96, and 192 hours, was used to determine 
the relative intensities given in Table II. The 
problem of selecting a crystal of the proper size 
to give reliable intensities is troublesome so that 
before using the experimental data it was neces- 
sary to examine the factors influencing the 
intensities in order to ascertain how exactly the 
sample used may be considered to have the 
properties of (1) a small nonabsorbing crystal, 
(2) a small absorbing crystal, or (3) a large 
absorbing crystal which intercepts the entire 
x-ray beam. For each of these three cases calcu- 
lation of intensity can be made with some relia- 
bility. The intensities given below are taken from 
photographs obtained with a crystal of dimen- 
sions approximately 0.10.1X0.03 mm. The 
assumption has been made that for this crystal 
absorption and extinction play a negligible role. 
This assumption is based on the following facts: 
(1) Lines of the same form reflecting at widely 
different angles of rotation of the crystal have 
the same intensity as nearly as can be deter- 
mined visually. (2) If the explanation of the 
observed intensity peculiarities is to be found 
in absorption or extinction effects no solution 
would result in many cases. This last argument 
is of a negative nature it is true, but it serves 
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to give validity to the principal conclusion to be 
drawn from this work, the tetragonal nature of 
the temperature motions of the silver atoms. The 
intensity comparison J 90) >J601) is an example 
of this last point and will be discussed later. 

The intensity calculations were made by 
means of the formula 


I~[(1+ cos? 26) /2 sin 26]| S\?7/? 


where |.S|? is the square of the structure factor 
and 7; is the temperature factor which will be 
discussed in detail later. The f values for silver, 
oxygen and phosphorus were taken from the 
tables of Pauling and Sherman.! 


Unit, SPACE GROUP, AND APPROXIMATE 
STRUCTURE 


The edge of the cubic unit has been found to be 
5.99;A to be compared with the value obtained 
by Wyckoff 5.99;A. The data are given in 
Table I. Laue exposures of 720 milliampere hours 
showed no reflections demanding a larger unit. 
Calculation from the unit size and the density 
gives the number of molecules in the unit as two. 
The observed absences (hkl, h=k, 1 odd) permit 
the space groups O,°, T,4*, O'!, O?, O,! and T 2. 
The last four may be ruled out by structural and 
intensity considerations as Wyckoff has shown. 
O,® may also be eliminated since it does not 
permit the separate existence of phosphate 
groups.* The special positions for T 4‘ are: 


2P at 000; 4, 4, 3. 
6Ag at 3, 3, 0; 5, 3, 0; 
80 at uuu, 

1 


2 
a+u, }- 


These conclusions are in good agreement with 
Wyckoff’s results. 

It will be noticed that the arrangement of the 
phosphate groups is body centered so that oxy- 
gen and phosphorus will not contribute to 
reflections with h+k+/ odd. In addition to the 


* The intensity inequality Z(so. >Jus1:) shows that the 
oxygen parameter must be larger than 0.30 thus eliminating 
the possibility of O,?. On a symmetrical Laue photograph 
where these two planes reflect at the same values of sin 0 
and ) the ratio of intensities is about four whereas for 
O,* (w=0.25) these reflections have the same structure 
factor. 

‘Linus Pauling and J. Sherman, Zeits. f. Krist. 81, 1 
(1932). 
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TABLE [. 








(hkl) sin 0 dhkl) (exp.) d(nkt) (cale.) 





2.997 
1.896 
1.498 
1.340 
1.176 
1.113 
0.985 
0.894 
0.823 
0.603 


2.989 
1.896 
1.498 
1.343 
1.179 
1.114 
0.985 
0.893 
0.823 
0.602 


(200) 
(301) 
(400) 
(402) 
(501) 
(502) 
(601) 
(603) 
(702) 
(804) 


0.1189 
.1870 
.2387 
.2640 
.3050 
3183 
.3680 
3792 
.4306 
5285 








regular space group absences the special posi- 
tions for the silver atoms lead to the following 
absences characteristic of this structure type: 
(hkl) with h, k, and / all odd; (hkl) with h—k 
=4n(n=1, 2, 3 etc.) and / odd. 

The absences and intensities can be accounted 
for in a rough qualitative way on the basis of the 
structure described above for “=0.35. Large 
discrepancies persist, however, which cannot be 
due to inaccuracies in the interpretation and 
which make it necessary to attempt some 
refinements. 


DISCREPANCIES AND MorRE ACCURATE 
> STRUCTURE 


Long Laue exposures with the x-ray beam 
making an angle of 7° with the axis [100] show 
the presence of a reflection from (401) very 
definitely. As shown above the appearance of 
this line is not compatible with the special 
positions of 7; assigned the silver atoms. It does 
not seem reasonable to attribute the appearance 
of the reflection to the oxygen atoms for there is 
no apparent physical reason for moving them 
from the three-fold axes or placing them in non- 
equivalent positions. Furthermore lines of the 
type (401) have been found on Laue and oscil- 
lation pictures at angles where the scattering 
power of oxygen is too small to account for the 
observed effects. 

Another type of disagreement between experi- 
mental intensities and those calculated for the 
above structure is shown by the comparison of 
the intensities of reflection from (600) and (601). 
The line (600) is nearly twice as strong as (601) 
and it reflects at about the same angle of crystal 
rotation so that absorption cannot be very 
important. The calculations show that (600) at- 
tains its maximum intensity at 1=0.33, since 
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here all the oxygen atoms reflect in phase, and 
that J 600)<J;601) even at this value of the 
parameter. It can be seen that neither absorp- 
tion nor an incorrect oxygen parameter can be 
responsible for this disagreement. 

A number of discrepancies such as this and the 
appearance of reflections from (401), (801), (621) 
and a doubtful very faint (841) show that for no 
value of the oxygen parameter can the calculated 
intensities even qualitatively be brought into 
agreement with those observed. It also seems 
clear that experimental conditions cannot be 
responsible for these effects. 

With the random. distribution assigned to 
silver in some previous work in mind, and be- 
lieving that there might be some tendency for the 
silver atoms to form only two bonds to oxygen 
(as in AgeO) I first attempted to get agreement 
by assigning two possible positions to each silver 
atom one on each side of the position 4, }, 0 along 
the twofold axis. To do this twelve equivalent 
positions of 74! were used and a new parameter, 
p, introduced which gave the separation between 
the new positions and the original one. The ap- 
pearance of lines of the type (401) could be ac- 
counted for in this way but an attempt to fit 
the intensities of the other lines failed. In order 
to force agreement with the experimental data it 
was found necessary to use eighteen positions, 
the original six and the twelve new ones and to 
weight these positions. The best agreement was 
obtained when the relative weights of the six 
and twelve positions were 2:1, p=0.035 and 
u=(0.34. In this case it may be added the agree- 
ment was only qualitative. 

This solution of the difficulty seems unsatis- 
factory and is, indeed very difficult to interpret 
physically ; it is included here only to show the 
analogy between this case and that of silver 
iodide. It was noted that the distribution along 
the twofold axes corresponding to this solution 
(with isotropic thermal motion about the points 
assigned to the fractional silver atoms) it is not 
greatly different from the probability distribu- 
tion of an anisotropic oscillator, the various 
states of which are weighted by an appropriate 
Boltzmann factor. This suggests that the silver 
atoms may not behave as spherically symmetrical 
oscillators but may have a larger amplitude along 
the twofold axes, i.e., in the directions of the new 
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positions which above have been given the com- 
bined weight one-half. 

The introduction of this sort of oscillator was 
accomplished by replacing the temperature 
factor ordinarily used, of the type 7; =e78in#? 
by a new factor derived on the assumption of an 
harmonic oscillator with tetragonal symmetry. 
By using the simple classical method of Darwin’ 
the expression for the temperature factor is 
found to be: 


Ty = eB (sin 8)? /d2 9B (B-1) (sin)? /d? cos? y 


= g—8 (sind)? /d2——B(B—1) (hi/ao)” 


where £ is the exponential coefficient which the 
factor would have if the force constant along 
the tetragonal axis were equal to the force 
constants perpendicular to this direction, B 
is the ratio of the force constant, perpendicular 
to the tetragonal axis to that parallel to it and 
cos y is the cosine of the angle which the fourfold 
axis of the oscillators makes with the normal to 
the reflecting plane. Simple substitution for sin @, 
\ and cos’ y gives the last expression which is in 
somewhat more convenient form. Here h; is the 
Miller index of the plane along the axis parallel 
to the tetragonal axis of the oscillator. In the 
case of this crystal, of course the formula must 
be applied to each set of two similarly oriented 
silver atoms in the unit. This is quite simple in 
practice. A table of e~8(8-)(*i/00)” against h; can 
be set up and the value of e~8(8-)) (*i/20)” can then 
be quickly multiplied by the f value obtained 
from the usual f vs. sin @ curve to obtain the 
effective scattering factor, since the f value from 
the curve may easily be computed containing the 
factor eB (sind)? /»?, 

The problem is now resolved into one of three 
parameters, u, 8, and B. By means of lines inde- 
pendent of u rough values for 8 and B were ob- 
tained. Lines of this type while showing some 
qualitative discrepancies for the case B=1 are 
not very sensitive to the actual value of B and 
more accurate values were obtained by con- 
sidering all lines by the trial and error method. 
The evaluation of the parameters was carried 
out by plotting J} against u for a series of values 
of 6 and B. In this way the following values 
were obtained: B=1.75+0.25, B=2.12+0.12, 
u=0.345+0.0025. 


5C. G. Darwin, Phil. Mag. 27, 315 (1914). 
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TABLE II. 








PLANE I(obs.) I(calc.) I(cale. 
u =0.345 
B=2. 
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* The intensities of these two lines were obtained from long exposure 
pictures over a small angular range so that the fact that they are 
calculated weaker than some lines not observed is not significant since 
these other lines were not in position to reflect on the long exposure 
pictures. 


The data are listed in Table II where the 
fourth column gives the best calculated intensi- 
ties for the case of spherically symmetrical 
oscillation of the silver atoms (B=1). The agree- 
ment will be seen to be well within the experi- 
mental error. In the case of some strong lines, 
(320), (400), and (404), the error is quite large 
but here the qualitative relationships are all 
correct. In the case of some of the weaker lines 
qualitative disagreements exist, but the difficul- 
ties of comparison are here so great that the 
discrepancies between experimental and calcu- 
lated intensity are again within the error of the 
determination. At large angles the lines become 
somewhat diffuse so that it becomes extremely 
difficult to make an accurate comparison of 
integrated intensities. It is true, too, that the 
total neglect of absorption is not justified. 
Considered as a whole, however, the agree- 
ment between experimental and calculated in- 
tensities is probably better than one would 
anticipate. 

The limits’ of error for 8 were established by 
making calculations for 8=1.5 and B=2.0 and 
observing that the best fit in these cases was 
markedly inferior to the agreement obtained 
with 6=1.75. The limits of error for B were 


fixed in the same way, calculations being made 
for B=2.00 and B=2.25. The lower limit for the 
oxygen parameter is quite definitely fixed by 
four reliable comparisons, the upper limit, 
however, is somewhat more uncertain. At 
u=0.3475 I 200) =T 400); the observed intensity 
relationship is I (400) =1.5 XJ 200). In this case it 
must be true that absorption has a considerable 
effect in reducing the relative intensity of (200). 
The upper limit for uw is assigned in the belief 
that the absorption for this small crystal cannot 
cut down the intensity by more than one-third. 


CONFIGURATION OF THE PHOSPHATE GROUP 


The special positions of T44 assigned to oxygen 
and phosphorus give the configuration of the 
phosphate group as a regular tetrahedron inde- 
pendent of the oxygen parameter. The parameter 
value (w=0.345+0.0025) leads to a phosphorus- 
to-oxygen distance of 1.61+0.03A, indicating 
that the P—O separation in this crystal is 
slightly larger than in KH2PQOx,, in which the 
potassium-to-oxygen bond is almost certainly 
ionic. The qualitative considerations given below 
will show that this effect is reasonable. The 
phosphorus-to-oxygen distance in the phosphate 
ion is influenced by some double bond formation 
and some ionic character of the P—O bond itself. 
Both of these effects presumably tend to de- 
crease the interatomic distance from 1.76A, the 
sum of the tetrahedral single bond radii accord- 
ing to Pauling and Huggins,® to the value 1.56A 
observed by West. Appreciable covalent bond 
formation between the oxygen atoms of the 
phosphate group and the surrounding silver 
atoms (see Fig. 1) will reduce the number of 
electrons on the oxygen atoms available for 
double bond formation and consequently the 
amount of double bond character of the P—O 
bonds should decrease. In accordance with the 
work of Pauling, Brockway and Beach’ this 
decrease in the amount of double bond character 
should show up in an increase in the interatomic 
distance, and furthermore, in a case such as that 
under consideration where here the bond has 


6L. Pauling and M. L. Huggins, Zeits. f. Krist. 87, 205 
(1934). ; 

7 Pauling, Brockway and Beach, J. Am. Chem. Soc. 57, 
2705 (1935). 
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fifty percent or more double bond character, the 
above work shows that a slight change in the 
amount of double bond character will have only 
a very small effect on the interatomic distance. 
One may deduce that the change in double bond 
character will be small from the fact that it 
depends on the covalent character of the Ag—O 
link which seems to be very nearly ionic since the 
interatomic distance here found is 2.34A+0.01A, 
or very nearly the sum of the ionic radii when 
corrected for coordination number. The sum of 
the covalent radii is by contrast 2.18A. 

It is difficult to estimate the influence of the 
ionic nature of the P—O bond. The formation of 
a Ag—O covalent bond should increase the ionic 
character of the P—O bond while at the same 
time the decrease in double bond formation 
should diminish it. On the basis of the experi- 
mental data one must conclude that the effect of 
diminished double bond formation outweighs 
that of increased ionic character, and, as a 
consequence of this, that the normal structure of 
the phosphate ion is essentially covalent. 


THE TEMPERATURE VIBRATIONS OF THE SILVER 
ATOMS 


The value obtained for B, the ratio of the force 
constants perpendicular and parallel to the four- 
fold axes of the silver oscillators gives an am- 
plitude of vibration along the tetragonal axis 1.46 
times as large as the amplitude perpendicular to 
this direction, if one makes use of the assump- 
tions involved in the calculation of the tempera- 
ture factor, a harmonic oscillator for which the 
amplitude is inversely proportional to the square 
root of the force constant. Te check this result a 
simple calculation has been made in which the 
silver-oxygen bonds have been treated as springs 
with no resistance to bending. With the dimen- 
sions taken from the crystal this calculation leads 
to a value of B equal to 2.75, or a ratio of the 
amplitudes of vibration equal to 1.65. The 
agreement here is probably well within the com- 
bined errors of calculation and experiment. The 
calculation of the temperature factor is un- 
questionably oversimplified. Certainly a con- 
sideration of the restoring force for bending of 
the bonds (if there is appreciable covalent 
character) would lead to closer agreement. 
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DISCUSSION OF SILVER IODIDE AND BROMIDE 


In a previous paper’ it has been shown that the 
x-ray spectra of AgI at temperatures from 100°C 
to — 180°C may be accounted for by assuming a 
random distribution of silver atoms among 
five positions inside the tetrahedron formed by 
the surrounding iodine atoms; one position in the 
center and four displaced toward the four faces 
of the tetrahedron. This structure for silver iodide 
may be changed in accordance with the proce- 
dure used in the case of silver phosphate. Here, 
however, the use of Hooke’s law forces leads to 
spherical symmetry for the atomic vibrations, 
and the anomalies of the x-ray spectra cannot be 
explained so simply. It is necessary to introduce 
fourth power terms into the potential function, 
and the calculation of the corresponding tem- 
perature factor becomes quite cemplicated. If 
we assume fourth power terms giving the poten- 
tial functions for the silver and iodine atoms 
tetrahedral symmetry we may argue by analogy 
with the silver phosphate crystal that these 
potential functions must have a smaller slope, in 
the directions of the faces of the surrounding 
tetrahedra since it was in these four directions 
that fractional parts of silver atoms were moved 
to get agreement with experiment. We are now 
permitted to treat iodine and silver atoms alike 
which seems more satisfactory than the as- 
sumption of special positions for silver alone (it 
is true, however, that we cannot expect the two 
atoms to show the same degree of anhar- 
monicity). If the silver and iodine tetrahedra 
were similarly oriented, one would expect no 
marked effect other than a decrease in scattering 
power, but in the wurtzite structure they are 
oppositely directed along the c-axis. This means 
that the anharmonic character of the vibrations 
of each atom contributes to the anomalies ob- 
served and that the differences in the amplitudes 
of vibration need not be very great to account for 
the observed effects. 

Photographs taken of AgI at low tempera- 
tures (— 180°C) indicated that the structure was 
essentially that of wurtzite with the ideal value 
of the parameter. Two reflections, (10-9) and 
(00-12) were in disagreement with this conclu- 
sion, (00-12) by a factor of ten. The assumption 
was made in the previous paper that of the 
five potential minima in each silver tetrahedron 
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the central one had the lowest energy so that 
with decreasing temperature the atoms would 
tend to settle in this position, but that at the 
temperature of the experiment only about four- 
fifths of them had done so. The explanation of 
the temperature effects becomes much simpler on 
the basis of the structure proposed here, and 
one would expect the discrepancies to persist 
even to 0°K because of the existence of zero- 
point vibrations, though they would naturally 
be small. 

The structure to be assigned to silver iodide on 
the basis of these considerations is that of ideal 
wurtzite in which the silver and iodine atoms 
behave as anharmonic oscillators having larger 
amplitudes in the directions of the faces of the 
tetrahedra formed by the nearest neighbors than 
in the directions of the corners of the tetrahedra. 
The fact that the effect is observed for this 
crystal and not for others with the same type 
structure is due probably to the fact that in this 
case of a very soft crystal the temperature factor 
and so the amplitudes of vibration are so large 
that the effects outlined above can be detected 
by the use of relative intensities obtained photo- 
graphically. 


An attempt has been made to find a similar 
effect for AgBr. That such an effect should exist 
seemed reasonable from the previous explanation 
of the structure of the iodide. No discrepancies 


between the observed intensities and thease 
calculated for the ideal sodium chloride structure 
have been found. The temperature factor, 
8=3.5 is rather high but what one would expect 
for such a soft crystal. If, as in the two previous 
cases we are to believe that observed discrepan- 
cies are to be attributed to deviations of the 
temperature motions from spherical symmetry, 
we can see quite easily that in the case of AgBr 
none should arise, for, not only would it be 
difficult to distinguish an oscillator with octa- 
hedral symmetry, the symmetry of the environ- 
ment of Ag and Br atoms (or ions), from one of 
spherical symmetry, but in this case also the 
octahedra have the same orientation with re- 
spect to the reflecting planes so that no large 
effect would be observed in any case. If this 
analysis is correct one must attribute the devia- 
tions from additivity of ionic radii in the inter- 
atomic distances of the silver halides, which are 
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intermediate between the sums of the covalent 


‘and ionic radii to an intermediate character of 


the bonds or, to what is equivalent, to polariza- 
tion. 

The phenomenon outlined above should be 
found wherever the environment of an atom is 
greatly different from spherically symmetrical. 
A clear cut observation of the effect in x-ray 
investigations will be dependent upon the 
orientation of the potential functions of the 
atoms with respect to the reflecting planes. This 
should be clear from the discussion of the cases 
of silver iodide and bromide. The effect will 
also vary in degree dependent on the mag- 
nitude of the force constants since it depends 
to a large extent on 8 as will be seen by refer- 
ence to the temperature factor developed 
for the case of silver phosphate. The aniso- 
tropic character of atomic vibrations should 
show up strikingly in soft crystals containing 
atoms with square coordination, or atoms which 
form two opposed bonds. PdO and KePtCl, are 
examples of the first class. For the first com- 
pound the data are insufficient for accurate 
analysis but in the second case some discrepan- 
cies between calculated and observed intensities 
in Dickinson’s work on K2PtCl,’ can be explained 
on the basis of tetragonal symmetry for the 
oscillations of the platinum atoms. The cuprite 
type structure is an example of the second class 
but here it is easy to see that the orientation of 
the O—Cu—O bonds makes it probable that the 
detection of the anisotropy of motion would be 
difficult to detect without accurate intensity 
measurements. It is probably true, too, that in 
many crystals where the configurations would 
allow the observance of the effects outlined the 
temperature oscillations are of such small 
amplitude at toom temperature that the de- 
parture from spherical symmetry would not have 
a great influence on the x-ray spectra. 

I am indebted to Professor Linus Pauling for 
discussion of this work and helpful suggestions. 

Note added in proof: Since the writing of this paper two 
articles treating the subject of anisotropic thermal motions 
of crystals have come to my attention, a theoretical dis- 
cussion of the phenomenon by Zener in the Physical Review 
of January and a note on the observation of the effects for 


the cases of Cd and Zn by G. W. Brindley in the February 
issue of Nature. 


®R. G. Dickinson J. Am. Chem. Soc. 44, 2404 (1922). 
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This section will accept reports of new work, provided these 


are terse and contain few figures, and especially few halftone. 


cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


The Conductance of Salt Crystals 


Recently Rodebush and Cooke! criticized a paper of the 
author? in a way which hardly seems adequate. 

(1) It is commonly known, as assumed in my article, 
that generally the conductance cannot be represented 
over the whole range of temperature by a single exponen- 
tial. It was the aim of my investigation to show that for 
regions, where a simple exponential holds—this is always 
the case for unipolar conductors and sufficiently high 
temperatures—one can derive the correct formula with 
constants of the right order of magnitude. Then it is not 
difficult to generalize the result for more complicated cases. 
Of course one can derive the same formal expression by 
more or less vague analogies with reaction rates as has 
been shown sufficiently often since the time of Rasch and 
Hinrichsen.* 

(2) Rodebush and Cooke’s attempt of a derivation of 
the self-diffusion formula is not without contradiction 
within itself. It holds generally that properties belonging 
to a stationary state—and self-diffusion as well as electro- 
lytic conduction refer to a stationary state—cannot be 
derived by analogies with the rate of a reaction without 
an explicit consideration of the reverse reaction. Rodebush 
and Cooke’s reasoning contains instead of the reverse 
reaction only the loss of mobility of an activated ion which 
moved away from its former lattice point. Therefore they 
describe the crystal in a state of increasing perturbation 
(‘‘Fehlordnung’’), not considering where the dislocated ions 
shall remain. A consequent treatment, however, would 
require as reverse reaction the return of an activated ion 
into a normal lattice position; therefore one first has to 
assume that vacant lattice points are present. If Rodebush 
and Cooke had consequently carried through their reason- 
ing they would have arrived at the conclusion that in 
equilibrium a certain fraction of the ions must be dislocated 
to irregular positions, leaving empty lattice points at their 
former positions; but this is just the model used by 
Frenkel‘ and the author. In special cases such empty 
positions are present in the original lattice structure as 
in the cases of a—Ag J® and a—Ag, S* which have been 
found after the author’s first paper. 

(3) Usually the methods of statistical mechanics and 
thermodynamics give more precise results in the treatment 
of states of equilibrium than in kinetic problems. Espe- 
cially for the treatment of disturbances of a crystal lattice 
(“Fehlordnungs-Erscheinungen”) the existence of which 
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can be concluded from numerous observations we refer 
to the fundamental papers of Schottky and Wagner.’ We 
mention, however, that the author’s assumption of the 
existence of an equivalent number of ‘“‘holes’’ and dislocated 
ions is not always true; for instance it may happen that 
the number of dislocated ions is negligible compared with 
the number of ‘‘holes’’ of both signs, which probably holds 
in the case of alkali halides.’:° 
W. Jost 
Institut f. physik. Chemie der tech. Hochschule, 


Hannover, Germany, 
March 25, 1936. 


1 W. H. Rodebush and T. G. Cooke, J. Chem. Phys. 3, 834 (1935). 

2 W. Jost, J. Chem. Phys. 1, 466 (1933). 

3 E. Raschand F. W. Hinrichsen, Zeits.f. Elektrochemie 14, 41 (1908). 

‘I. Frenkel, Zeits. f. Physik 35, 652 (1926). 

5.L. W. Strock, Zeits. f. physik. Chemie B25, 441 (1934). 

6 P. Rahlfs, Zeits. f. physik. Chemie B31, 157 (1936). 

7C. Wagner and W. Schottky, Zeits. f. physik. Chemie B11, 163 
(1930); C. Wagner, Zeits. f. physik. Chemie, Bodenstein Festschrift, 
(1931), p. 177; Zeits. f. physik. Chemie B22, 181 (1933); C. Wagner, 
Physik. Zeits. 36, 721 (1935). 

8 W. Schottky, Zeits. f. physik. Chemie B29, 335 (1935). 

9cf. W. Jost, Zeits. f. tech. Physik 16, 363 (1935). 





The Raman Spectrum of Oxalic Acid 


In a recent note Angus and Leckie! have taken exception 
to an observation of Hibben made in a previous publica- 
tion.2 This article recorded and discussed the Raman 
spectrum of oxalic acid in aqueous and alcoholic solutions, 
as the dihydrate and as the anhydrous compound. It was 
observed, in comparison with other organic acids, that the 
two C=O Raman shifts (Av 1645 and 1740) obtained from 
anhydrous acid and aqueous solution are anomalous. In 
alcoholic solutions the doubling of the carbonyl shift no 
longer occurred and one shift, Ay 1755, resulted, which is 
more in accord with the carbonyl shift obtained from an 
ester, aldehyde or ketone. Consequently it was affirmed 
that, based on the magnitude of the two frequency dis- 
placements, the two carboxyl groups behave differently in 
aqueous solution and in the anhydrous acid, but iden- 
tically in alcohol. The attribution of Ai 1430 to an O—H 
shift in oxalic acid by Rao,* Tschitschibabin’s and the 
hexahydroxy formulas were discussed in detail. 

In addition to the above, which comprised the greater 
portion of the article in question, it was stated: ‘‘No shifts 
corresponding to those at approximately Av 1650 and 1750 
of the anhydrous acid were observed in the hydrate. This 
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does not preclude the possibility that they may be extant, 
but if so they are very much weaker than the same shifts 
in the dehydrated acid.” It is principally with this remark 
that Angus and Leckie disagree, as they state: ‘‘. . . the 
Raman spectra obtained from the saturated aqueous 
solutions of oxalic acid and from the dihydrated oxalic 
acid crystals are substantially the same.’’ However, they 
add: ‘‘Our results, obtained from crystals which were 
definitely dihydrate crystals, agree with those of Rao . . .” 
and, ‘‘Since Hibben’s conclusions that the two carboxyl 
groups behave differently depend to a large extent on the 
difference in the spectra of the anhydrous and hydrated 
crystals we submit that such conclusions may be in- 
validated by his comparing the spectrum he obtained for 
his anhydrous crystals with Rao’s spectrum for the solid.” 

The last quotation leaves such an erroneous impression 
that it is necessary to correct it. As pointed out in the 
first paragraph of this note, the different behavior of the 
two carboxyl groups is predicated on the observed magni- 
tudes of the two carbonyl shifts and not on a comparison 
of the anhydrous, or hydrated material with the solution, 
consequently these conclusions are not affected by the 
presence or absence of such shifts in the hydrated material. 

The apparent weakening of these lines in the anhydrous 
acid and their greater weakening (or absence) in the 
hydrated compound were explained by this writer as 
possibly attributable to chelation. If Angus and Leckie’s 
observations are correct this explanation would be without 
significance. Nevertheless Rao’s results cannot be cited 
in entire support of themselves since his published photo- 
graphs show clearly that in the hydrate or anhydrous acid 
(it is immaterial to the argument which he used) the 
relative intensities of the carbonyl frequencies are markedly 
reduced as compared with the solution, in contradistinction 
to the first statement of Angus and Leckie’s previously 
quoted. The intensities published by Rao are: (for A> 1640 
and 1740, respectively) crystal: 00, 1; solution: 1, 3; or 
approximately a 66 percent decrease taking Ay 845 as a 
reference standard. It is unfortunate that in republishing 
Rao’s data Angus and Leckie omitted the intensities, 
which indicate an agreement with the present writer. 

The fact that carbony] shifts may weaken is already well 
known from the spectra of acid salts. So far as the hydrate 
is concerned Zachariasen* from x-ray measurements has 
concluded that the C—O bond in this compound is at 
least part of the time a single bond, and that this substance 
is not a normal hydrate of oxalic acid. 

It is helieved that some of the difficulties of interpreta- 
tion may have arisen in the use of an unfiltered mercury 
spectrum as a source of excitation. The O—H band, un- 
usually sharp, falls in approximately the same position as 
the C=O shifts and the 4730 Hg line overlaps in the same 
position. 

James H. Hippen 

Geophysical Laboratory, 


Carnegie Institution of Washington, 
January 21, 1936. 


1 Angus and Leckie, J. Chem. Phys. 4, 83 (1936). 
2 Hibben, J. Chem. Phys. 3, 675 (1935). 

3S. Rao, Zeits. f. Physik 94, 536 (1935). 

* Zachariasen, Phys. Rev. 45, 755 (1934). 


THE EDITOR 


Reply to Dr. Hibben’s Letter 


Dr. Hibben’s criticisms may be answered very briefly 
His suggestion that our results for the dihydrate crystals 
are vitiated by incomplete filtration of the 4047A Hg line 
is incorrect since we used a strong solution of NaNO: as a 
filter and the 4047 undisplaced Hg lines were very weak 
on the plates. His second quotation from our letter tends 


‘to misrepresent the arguments put forward by us. We 


must maintain our view that Rao’s crystals, showing 
frequencies at approximately 1660 and 1760 cm™', were 
dihydrate crystals and that Rao’s results, therefore, sup- 
port our own observations. We have satisfied ourselves 
that our results are correct. 

The purpose of our letter was merely to indicate dif- 
ferences in results which, as we pointed out, may invalidate 
some of Hibben’s conclusions. We do not think that the 
point at issue is sufficiently important to merit an extended 
controversy. 

W. RoGiE ANGUS 
ALAN H. LECKIE 
University College, 


London, W.C.1., 
February 19, 1936. 





Metallic Binding According to the Combined Approxima- 
tion Procedure 


A combined»approximation procedure has been pro- 
posed by one of us! for the calculation of the chemical 
interaction between atoms. In this method the atomic 
kernel is treated according to Thomas-Fermi and the 
valence electron by a Schrédinger equation. The total 
action of the kernel upon the valence electron—inclusive 
of the Pauli principle—then goes over into an additional 
potential in the Schrédinger equation, which can be taken 
approximately from the spectra. The application to Kk; 
and KH gave plausible results, however, on account of 
the omission of higher approximations they were still not 
absolutely convincing. We have therefore applied the 
combined approximation procedure to the problem of 
metallic binding. Here one can easily take into account 
the interactions between the valence electrons to a higher 
degree of approximation than in single molecules. 

It was established in the beginning that the law 
—1/r+(A/r)e*" for the total field of the kernel gives a 
good approximation in the case of the free Mg atom: When 
A and « were so determined that the lowest s term and 
p term of Mg* were rendered correctly, then Ritz's method 
(with riz in the eigenfunction) yielded good results for the 
ground term of Mg. For the alkali metals A and « were 
determined from the two lowest s terms with the aid of 
one-parameter hydrogen-like eigenfunctions. Since the 
field of the ion determined in this way runs a rather flat 
course, it is sufficient in the first approximation to set the 
valence electron in the lattice equal to a plane wave. The 
interaction of the kernels arising from their overlapping 's 
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neglected. The energy of the metallic lattice is then com- 
posed of the following parts: 


(1) Electrostatic energy ~—1/d (taken from Wigner and 
Seitz?) 

(2) Exchange energy after Bloch? ~—1/d 

(3) Correlation energy after Wigner‘ (from Fig. 7 of 4) 

(4) Thomas-Fermi energy ~1/d? 


a , , 1 
(5) Additional energy = El. density x f e 
r 


The Table I gives the values found in this way for the 
heat of sublimation and do*. Through improvement of the 
uniform distribution of the valence electron in the lattice, 
a further approach of the heat of sublimation to the experi- 
mental value may be attained. That the additional energy 
reproduces well the influence of the kernels upon the 
valence electron is demonstrated by the calculated @?E/aV? 
which agrees much better with the experimental data 
than the calculation of Rice,> who instead of our intro- 
duction of the additional energy assumed the kernels to be 
impenetrable and determined their radius empirically 
from do itself. The relationship between sublimation energy, 
ionization energy and the work function first established 
by Wigner and Bardeen® is also valid in this theory. 


TABLE I. 








@E/dV2 (E: EN- 
ERGY, V: ATom- 
1c VOLUME) 
(1083 erg-cm~) 


Work Func- 
TION (ev) 


BINDING 
ENERGY (ev) 


LATTICE 
Const. (A) 





calc. 


cale. exp. calc. exp. calc. exp.5 Rice5 calc.  exp.§ 





4.35 4.24 
5.20 5.25 
5.75 5.62 
6.15 6.05 


0.91 1.13 1. 1.96 2.25-2.46 
0.85 0.94 0. 1.95 2.17-2.24 
0.594 0.865 0. 5 / 1.90 2.16-2.19 
0.599 0.83 0. 1.86 1.81-1.96 








Column 4 gives the work function calculated in this way, 
omitting the double layer term. When one bears in mind 
that the only empirical data used are the spectra of the 
free atoms the agreement is satisfactory and there is the 
hope that this simple approximation method will prove 
useful for further metallic problems where the more 
vigorous methods become involved in formal difficulties. 
The complete work will appear in the Acta Physico- 
chimica U.R.S.S. - 
H. HELLMANN 
W. KAssATOTSCHKIN 
Karpow-Institute for Physical Chemistry, 


Moscow, 
March 15, 1936. 


*Li was omitted, as in this case the application of a statistical 
method is not free of arbitrariness. 
_1H. Hellmann, J. Chem. Phys. 3, 61 (1935); Acta Physicochimica 
U.R.S.S. 1, 913 (1935); 4, 225 (1936). 

? E. Wigner and F. Seitz, Phys. Rev. 46, 509 (1934). 

°F. Bloch, Zeits. f. Physik 57, 545 (1929). 

4E, W igner, Phys. Rev. 46, 1002 (1934). 

50, K. Rice, J. Chem. Phys. 1, 649 (1933). 

‘E. Wigner and J. Bardeen, Phys. Rev. 48, 84 (1935). 
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The Raman Spectra of Carbon and Silicon Tetrafiuorides 


We have photographed the Raman spectra for the tetra- 
fluorides of carbon and silicon in both the liquid and 
gaseous states. The results are presented in Table I 
together with the Raman frequencies of CCl, for purposes 
of comparison. The relative intensities are given in 
parentheses. 


TABLE I. 

we (cm~!) w3 (cm~!) 
437 (1) (1200) 

( 


285) (1000) 
214 775 


ws (cm*!) 
635 (1) 
(431) 
311 


w1 (cm™~!) 
904 (10) 
800 


450 


CF4 
SiF« 
CCh 

So far only w; for Sif, has been found, the other three 
lines being apparently of small intensity. For gaseous CF, 
only w was definitely seen on the plates, and it was sharp 
and strong. w: and w, together with w, were obtained from 
CF, (1) at —120°. 

The CF, was prepared by treating carbon (Norit) with 
fluorine gas and submitting the raw product to a number 
of fractional distillations. It is difficult to prepare pure CF,. 
The SiF, was prepared by a conventional method. 

The assignment of frequencies is based on a number of 
considerations. Since w is the completely symmetrical 
vibration, it should appear intense and sharp in the Raman 
spectrum. Accordingly the strongest and sharpest ob- 
served line was assigned to w. It was found, from a con- 
sideration of the Raman spectra for eight tetrahalides 
that certain regularities existed (approximate constancy of 
ratios of corresponding frequencies), and on the assumption 
that no discontinuities in these regularities would occur 
with CF, the frequencies w: and w,; were assigned. 

The frequencies given in parenthéses were computed 
from empirically found ratios, or from the formulae derived 
from normal coordinate treatments of tetrahedral mole- 
cules by Urey and Bradley,' and Rosenthal.? 

Two absorption bands in the infrared spectrum of CF, 
have been observed recently by Eucken and Bertram,’ 
namely, w3;=1350 cm™ and w;=653 cm. The latter is in 
fair agreement with 635 cm~ found in the Raman spec- 
trum, and the estimate (1200) for w; is, considering the 
theory involved, in accord with w;=1350 cm™. Eucken 
and Bertram estimate w: to be 503 cm™ from specific heat 
measurements, and this agrees only moderately well with 
the observed value, w2.=437 cm. The disagreement is 
probably to be assigned to errors in the theory or measure- 
ment of the specific heats. 

Don M. Yost 
Epwin N. LASSETTRE 


SIEGFRIED T. GROSS 
California Institute of Technology, 
Pasadena, California, 
April 6, 1936. 


1 Urey and Bradley, Phys. Rev. 38, 1970 (1932). 
? Rosenthal, Phys. Rev. 45, 538 (1931). 
3 Eucken and Bertram, Zeits. f. physik. Chemie B31, 361 (1936). 





